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enhance the overall power density and efficiency. Tnis 
thesis is concerned specifically with the design and testing 
of the electrodynamic transducer which 1S responsible for 
the electro-acoustic power conversion. A computer model or 


e 


the driver/resonator system is presented along with the 
techniques for measurement of the electrical and mechanical 
parameters used as input for the model. 

A final driver design (including dimensional drawings) 
mea iilzing a modified EL 2450J neodymium-iron-boron 
compression driver and associated leak-tight electrical 


feed-throughs, microphone, accelerometer, pressure gaze, 


pressure housing, and resonator interface is provided. 
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A. BACKGROUND 
IE HIStoODy 

Thermoacoustics can generally be described as tne 
study of the interaction between heat and sound. Scientific 
interest in this interaction is not new. Lord Rayleizh (1878 
and 1945) discussed various qualitative examples of the 
production of sound oy heat. In one of these examples he 
describes investigations by Sondhauss (1850) of an effect 
observed by giassblowers: heating a bulb of glass at the 
end of a hollow rod occasionally led to the generation of 
sound. Faraday (1818) demonstrated that this effect occurred 
With gases other than air. Sondhauss proved that the vibra- 
tion of the glass itself did not generate the sound, but he 
offered no explanation as to what did. In his description of 
Sonhauss' work, Rayleigh stressed the importance of the 
Phase difference between temperature and particle motion. 
Rijke (1859) built an open tube with a wire mesh inside. 
When the mesh was heated, the tube produced a sound of 
impressive intensity. The functioning of this oscillator is 
slightly different since it requires steady gas flow for its 
operation. A similar sound production effect was observed by 


Taconis, et al. (1949) in hollow tubes immersed in liquid 


1 / D y ; 1 3 
helium (4.297) The Taconis oscillations were addressed 


quantitatively by Yazaki, et al. (1980). The Work descrites” 


above deals primarily with prime movers--devices that con- 
vert a temperature gradient to sound energy. We call this 
the classical era of thermoacoustics. 

Interest in thermoacoustics was renewed when the 
idea of the use of acoustical devices as refrigerators 
appeared. Gifford and Longsworth (1966) described the pump- 
ing of heat along a Surface caused by a periodic change in 
the pressure of the adjoining gas. Such a change can be 
produced by the oscillations of a sound wave. In their 
experiments Merkli and Thomann (1974) explored the heating 
and cooling effects on the wall of a gas-filled resonant 
tube. They found that heat was transported from a region 
near the velocity antinode (or maximum value) of the sound's 
standing wave to regions near the adjacent pressure  anti- 
nodes. Thermoacoustic prime movers have also received recent 
attention. Kempton (1976) investigated the excess noise of 
aeroengines above that predicted by theory. He determined 
that the sound was produced by unsteady heat transfer. Each 
of these latter three groups of exprimenters used some 
theory for comparison, but it was mostly qualitative. None 
provided the complete quantitative theory that would explain 
their experimental observations. 

The theoretical breakthrough in thermoacoustics 
field was made by Nicklaus Rott (1969, 19 449 1975, and 


1980). His theory combined basic principles from physics, 


EJ 


thermodynamics, and acoustics to quantitatively describe the 
effects found in both types of thermoacoustic devices: prime 
movers and heat pumps. Prime movers use a temperature gra- 
dient to create sound, as in the effects discussed by Sond- 
hauss, Rayleigh, Taconis, and Kempton. Heat pumps, or refri- 
gerators, use the oscillating pressure of a sound wave to 
produce a temperature difference, as described by Gifford 
and Longsworth and Merkli and Thomann. Rott described the 
effect found by Merkli and Thomann as thermoacoustic stream- 
ing. It is this effect that makes thermoacoustic refrigera- 
tion possible. 

Inspired by Ceperley's (1979) traveling wave Ster- 
ling Cycle heat engine and Rott's quantitative theory, the 
team of Wheatley, Hofler, Swift, Migliori, and Garrett 
(1982, 1983a, 1983b, 1985, and 1986) developed a series of 
thermoacoustic experiments at Los Alamos National Laboratory 
in New Mexico. They investigated the basic  thermoacoustic 
effects in both prime movers and refrigerators and compared 
their experimental results to Rott's theory. 

2. Thermodynamics 

We'll digress from history here to explain the ther- 
modynamic distinction between prime movers and refrigera- 
tors. This discussion follows Sears and Salinger (1975). A 
Exe mover receives heat at a high temperature, T,, (the hot 
reservoir), does work on its surroundings, and rejects heat 


at a lower temperature, Tc (the cold reservoir), as shoun in 


W 


Figure I-1.a. The first law of thermodynamics tells us that 
the heat flow (Qu) from the hot reservoir must equal the 
work done on the surroundings (W) plus the heat flow (Qc) to 
the cold reservoir: 
Qu = W + Qe. (I-1) 

The second law of thermodynamics states that in every  pro- 
cess the entropy of the universe must either remain constant 
or increase. The entropy change for the isothermal heat 
transport to/from the reservoirs is defined as the heat flow 
divded by the reservoir temperature. For a prime mover the 
second law can be written: 

AS4 E Ass cr Gaya eee gees (iam 
The thermal efficiency of the prime mover is defined as the 
ratio of the work output to the heat input: 


mn = W/Qn = (Qu - 0Oc)/Om = 1 - Qc /Qu. (I-3) 


Comparing equations 1-2 and 1-3 gives us a limit on the 


efficiency based on the temperatures of the reservoirs: 


. 


nail = Tey ia (I-4) 

This limit is known as the Carnot efficiency. 
In a refrigerator the process is basically reversed 
(see Figure I-1.b). Heat (Qe) is removed from the low tem- 
perature reservoir, work (W) is received from the surroun- 
dings, and heat (Qu) is delivered to the higher temperature 
reservoir. The expression for the first law is the same as 
for the prime mover, Equation I-1. For the refrigerator it 


means that the heat flow into the hot reservoir must equal 


FIRST LAW: Op" W + O. 


Er 


SECOND LAW: 
Te 


ap 


COEFFICIENT OF PERFORAANCE: 
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The Thermodynamics of (a) Prime Movers 


Figure I-1. 
and (b) Refrigerators. 


the heat flow from the cold reservoir plus the work done to 
the system. The second law for this case is: 


AS. E ASc al Qu/ Tu E Qc/Tc. ( 1-99 


The efficiency parameter for the refrigerator is called the 
coefficient of performance (COP). It is defined as the ratio 
of the heat flow from the cold reservoir to the work input 
to the refrigerator, or: 


COP 2» Qc/W s Qc/(Qu-Qc) ¿2 Te/(Tum-Tc). (I-6) 


The limit is known as the Carnot coefficient of performance. 

The work in this thesis deals solely with .etrigera- 
EC Even though prime movers are mentioned from time to 
time, the primary focus of the rest of our discussion will 
be refrigerators. 

3. Overall Efficiency 

According to Hofler (1986), the early thermoacoustic 
work at Los Alamos focused on experimental refrigerating 
engines. The performance of these engines fell short of 
expectations, leading to simple experiments on basic thermo- 
acoustic effects and a proof-of-principle experiment on 
thermoacoustic refrigeration. Hofler then applied the Rott 
theory to the experimental systems and solved the resulting 
equations numerically. For his doctoral dissertation trom 
the University of California, San Diego, Hotler designed and 
constructed a completely functional thermoacoustic refriger- 
ator. He also made accurate measurements of its thermody- 


namic efficiency, and used this efficiency to make compari- 


sons to the Rott theory. After receiving his doctorate, 
Hofler came to the Naval Postgraduate School (NPS) as a 
post-doctoral fellow, and Drought his prototype refrigerator 
with him. 

The purpose of this thesis, imeem uMet Lon with the 
work of several other students, is to modify Hofler's re- 
frigzerator design in order to improve its overall efficien- 
cy and make it suitable for space cryocooler applications. 
The basic desizn of the thermoacoustic refrigerator is shown 
schematically in Figure I-2. The driver (A), which produces 
the sound waves, is coupled via bolts to the resonator (D) 
Via a reducer cone (B) and bellows (C). Inside the resonator 
is a stack of plastic plates (E) and their heat exchangers 
(F) which allow the heat to be removed from the hot end arid 
absorbed by the cold end. It is the interaction between the 
sound waves and the plastic plates that produces a tempera- 
ture difference across the plates and/or pumps heat from the 
cold heat exchanger to the hot heat exchanger. A brief 
description of efficiency would be useful here. Figure I-3 
shows a diagram of the chain of energy transformations which 
occurs in a thermoacoustic refrigerator. There are three 
conjugate variable pairs (V and I; P and U; and AT and Ó-). 
and two transformations (electrical-to-acoustical and acous- 
tical-to-thermal). 

We start with electrical input power to the driver, 


which can be calculated (using r.m.s. values) with: 


Figure 1-2. Basic Thermoacoustic Refrigerator Showing 
(A) Driver, (B) Reducer Cone, (C) Bellows, 
(D) Resonator, (E) Stack, and (F) Heat 
Exchangers. 
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Figure I-3. Chain of Energy Transformations in a 
Thermoacoustic Refrigerator. 


Parac - I V cos 9, (I-7) 
where I is the current into the electroacoustic driver, V is 
the voltage across the driver terminals, Əə is the phase 
angle between I and V, and cos ə is known as the "power 
factor." An electrodynamic driver converts the electrical 
power to acoustic power (the first transformation).  Analo- 
gous with the previous definition of the electrical input 


power, the acoustic power is given by: 


Ha P U Cos Wi (1-8) 
where P is the acoustic pressure, U is the volumetric velo- 


city which equals the product of particle velocity Wage 
resonator cross-sectional area, and * is the phase angle 
between them. The efficiency of this electroacoustic trans- 
formation 1S -e2iVen Dy: 
nea = W/Poiec. (I-9) 
In the second transformation, the refrigerator con- 
verts the acoustic power to a temperature gradient (7T) and 
heat flow (Q). As discussed previously, the efficiency for 
this transition is given as the coefficient of performance: 
COP = Gc/W (I=-TOR 
The overall thermoacoustic efficiency is therefore the pro- 
duct of the efficiencies for the two transformations: 


Ntra = Mea * COP = Qr Pae ec. (I-11 


Improvements in the overall refrigerator efficiency 
can thus be divided into two distinct, though interrelated, 


efforts. The first is the subject of this thesis: modifica- 


10 


E3oOns Pome thie dr Vento optimize the electroacoustic 
efficiency. The second effort involves the optimization of 
the coefficient of performance and is addressed by LT 


Michael P. Susalla, USN, in his master's thesis. 


B. SPACE THERMOACOUSTIC REFRIGERATOR (STAR) 
1. Motivation 

In addition to improving overall efficiency, our 
experimental refrigerator will be designed and built with 
the intention of flying it on the Space Shuttle as part of 
the National Aeronautics and Space Administration's (NASA) 
Get Away Special (GAS) program. As of 3 February 1988, there 
is a Memorandum of Agreement between the Naval Postgraduate 
School (NPS) and the Air Force which provides funding for 
the space flight and which assigns the NASA payload. number 
G-337 to this project. The reason for testing the thermo- 
acoustic refrigerator in space is that the immediate appli- 
cations of the STAR are expected to be the cooling of elec- 
tronics, high Te superconductors, and infrared detectors in 
Space. 

There are currently two primary methods for cooling 
satellite systems in space: (1) evaporation of expendable 
cryogens (liquid helium, nitrogen, ammonia-methane, and 
SOlid hydrogen. etc.), and (2) closed cycle refrigerators 
fotariine cycle, Vuilleumier cycle,  etc.), which involve 


large reciprocating piston masses operating at low 


p 


frequency. (Walker (1983) and Smith, et al. (1984)] The 
disadvantages of these current cooling methods are their 
short lifetimes (expendable cryogens) and high vibration 
levels and low reliability (closed cycle refrigerators). The 
advantages of the thermoacoustic refrigerator over these and 
other systems lies in its simplicity and reliability. The 
STAR has no sliding seals, efficient performance, low vibra- 
tion levels, and (hopefully) a longer life span. 

The thermoacoustic refrigerator needs to be surroun- 


ded by a vacuum and insulation material to minimize heat 
loss due to thermal conduction and radiation. The vacuum of 
space will provide the best insulation to determine the 
refrigerator's ultimate efficiency. The absence of gravity 
removes the possibility of thermal convection as an  addi- 
tional nuisance heat transport mechanism. 
2. Get Away Special (GAS) Program 

NASA's GAS program allows small, self-contained 
payloads to fly on the Space Shuttle in GAS canisters for 
relatively low cost ($10,000) (Get Away Special Team, 1984). 
The GAS can is five cubic feet in volume and can house a 
payload of up to 200 pounds. Each payload must contain its 
own electrical power, control, data acquisition and storage 
facilities, etc. The Shuttle'*s astronauts will turn 2 SN 
on or off at designated times during the flight, but are 


otherwise not involved with the experiment. 


B 


3. Shared Subsystems from NASA G-313 


Another group of NPS students and faculty (Boyd, 
ec al., 1987) have taken advantage of the GAS program to 
measure the resonant acoustic modes of the shuttle payload 
bay and the ambient acoustic environment produced as a 
result of main engine and booster operation during launch. 
This experiment is titled "The Space Shuttle Cargo Bay 
Vibroacoustics Experiment" and is designated by NASA as 
payload G-313. Several subsystems that were developed for 
NASA G-313 will be used by NASA G-337 (STAR). A schematic of 
the STAR in its GAS can is shown in Figure I-4. One of the 
Systems borrowed from NASA G-313 is the electronics system 
used to run the experiment and record the data. The recorder 
System consists of the INTEL model BPK SV75 magnetic bubble 
memory module and an NSC 800 microprocessor-based control- 
ler. Two other NPS students, LT Charles B. Cameron, USN and 
CPT Ronald Byrnes, USA, will be designing the analag elec- 
tronics and software to integrate these systems into the 
STAR experiment for their master's thesis. A block diagram 
of the electronics is shown in Figure I-5. 

Another borrowed system is the power supply. which 
consists of Gates brand lead acid battery cells (five 
ampere-hour, two volts each). These gelled electrolyte bat- 
teries are ideal for the STAR due to their high power den- 
ancy , low cost, and the absence of outgassing during dis- 


Charge cycles. NASA G-313 used a one layer battery of 68 
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Figure I-4. Schematic of the Space Thermoacoustic 
Refrigerator in its Get Away Special 


Cannister. 
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Figure I-5. Block Diagram of the Electronics for the 
Space Thermoacoustic Refrigerator. 
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cells providing 680 watt-hours of energy and weighing about 
80 pounds (including the cells’ support structure). Ne 
be using two battery layers with as many cells as we can, up 
to 136. This will give us a minimum of 680 watt-hours and a 
maximum of 1,360 watt-hours of available electrical energy, 
and a battery weight between 80 and 160 pounds. The number 
of battery cells we can use will depend on the total weight 
of the driver-resonator assembly and itS auxiliary equipment 
(vacuum can, gas reservoir, etc.) and electronics. 

The use of the GAS can imposes certain restrictions. 
Since we are using batteries to supply the power, the refri- 
gerator has to be energy efficient. Also, the GAS can setup 
requires the STAR to be compact and lightweight. Figure I-6 
shows a photograph of Hofler's prototype refrigerator. This 
setup is approximately six feet high. In comparison, the 
maximum payload height for the GAS can is 28.25 inches, or 
less than 2.5 feet. These considerations played a major role 


in the choice of equipment for and the design of the STAR. 


Qa. ¡SCORE 

Chapter II discusses basic acoustical and thermodynamic 
theory as applied to STAR. The results of this theory deter- 
mine the driver requirements. The theory of the selected 
electrodynamic driver is given next. Following the driver 
theory is a presentation and discussion of a computer 


program that models the coupled system, which consists of 
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Photograph of Hofler's (1986) 


Prototype Refrigerator. 


Figure I-6. 
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the  electrodynamic driver and the resonator containing the 
gas mixture and stack. 


Chapter III describes the technique for measuring the 


electrodynamic driver parameters. These parameters are used 
as input for the driver design, the computer model, and the 
efficiency calculations. This measurement technique was 


developed using an Altec 290-16K electrodynamic driver and 
was applied to the final driver. 


Chapter IV describes the measurement of the STAR driver 


parameters. The driver chosen for STAR is a custom unit, 
built by Harmon-JBL!‘™ _. which uses Mmeodimium-iron= ona 
(NIB) magnets. The STAR driver is a Minomemedift lication eum 


new line of compression drivers usinz NIB magnets and has 
characteristics similar toWNthe JBL 24599 

Chapter V discusses the construction of the STAR driver, 
including the housing and accessories (reducer, stiffener, 
bellows, accelerometer, microphone, electrical feed- 
throughs, etc.) necessary to integrate it with the resonator 
and the GAS cannister system. Dimensional shop drawings for 
the parts are included in Appendix C. 

Chapter VI gives conclusions and recommendations for 


further development. 
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EX LHERMOACOUSTIC THEORY 

Thermoacoustic theory has been developed in detail by 
Rott (1969, 1974, 1975 and 1980) and adapted to the thermo- 
acoustic refrigerator by Wheatley et al. (1982, 1983a. 
1983b, 1985 and 1986), Wheatley and Cox (1985), Hofler 
(1986) and Swift (1989). 

The space thermoacoustic refrigerator (STAR) basically 
consists of an acoustic driver producing sound waves in a 
resonant tube (see Figure I-2). This resonator is filled 
with a mixture of helium-xenon gas (12.5% xenon) pressurized 
to ten atmospheres, and contains a stack of plastic plates. 
This Chapter will present a qualitative model for the ther- 
moacoustic heat pumping process followed by a quantitative 
development for heat and work flow at a plate. The Chapter 
concludes with the efficiency of a stack of plates that are 
very much shorter than one quarter of a wavelength. 

1. Qualitative Picture 


Consider a parcel of gas that moves back and forth 


along one of the plates at the acoustic frequency (see 
Figure  II-1). As it moves, the gas parcel will experience 
changes in temperature and volume. Part of the temperature 


changes come from the adiabatic compression and expansion of 


the gas by the sound pressure, and part as a consequence of 
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Gas Parcel Diagram. (Note: Displacement 
in the Vertical Direction is Shown for 
Clarity--Parcel Moves in Horizontal 
Direction Only 


the local temperature of the plate itself. A temperature 
gradient may develop along the plate as a result of the 
operation of the refrigerator. The temperature and volume 
changes can be described by six separate steps (the changes 
in volume are indicated by the change in size of the square 
parcel in the figure). 

Assume the plate 1s at a uniform temperature, T. The 
parcel starts at position 1 (X = 0) with a temperature of T. 
The sound wave moves the parcel to the right to position 2. 
The parcel has undergone adiabatic compression and its tem- 
perature is now T**. Since the temperature of the gas is 


Bucher than that of the plate, heat (Q) will flow from th 


M 


gas parcel to the plate. The parcel's volume is decreased 


Ye 


(D 


ema its temperature lowered to T°. This heat flow alse 
causes the plate's temperature to increase at the position 
X*. Parcel position 3 is actually in the same location as 
position 2, but it is displaced vertically in the figure for 
Clarity. The parcel now moves to position 4 and expands 
EXcbatically to a new temperature of T- A repeat of this 
process puts the parcel in position 5 with a temperature of 


T777. Now the temperature of the gas is lower than that of 


the plate and heat (Q') flows from the plate to the gas 


parcel, expanding the parcel and raising its temperature to 
Ec This flow causes the plate's temperature to decrease at 
position xX-. The parcel moves to the right again under 


adiabatic compression and we are back where we started in 


position 1 with temperature T (adapted from Wheatley, et 
al., 1985). 

This gas parcel cycle is repeated at the resonator's 
operating frequency. For the STAR this frequency is between 
approximately 250 and 600 cycles per second, depending on 
the gas mixture (specifically the speed of sound in the gas 
mixture) and the length and shape of the resonator. 

Notice that the temperature of the gas parcel at 
position X = O is different depending on which direction the 
gas parcel is moving. It ıs this phase shift in temperature 
relative to motion that produces the thermoacoustic effect, 
as we will show in the calculations that follow. 

It is also important to note that in order tor mhen 
to flow between the gas parcel and the plate, the parcel 
must be vertically located within about a thermal penetra- 
tion depth of the plate. The thermal penetration depth (5) 


is the distance that heat can diffuse through the fluid 


during a time  1/«sm, where w is the acoustic angular 
frequency. 
Wheatley et al. (1986) describe thermoacoustic en- 


Eines as consisting of long trains of these gas parcels, all 
about a thermal penetration depth from the plate. The par- 
céls draw heat from the plate at one extreme of their oscil- 
latory motion and deposit heat at the other extreme. Adja- 


cent heat flows cancel except at the ends. The net result is 


bh 
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that an amount of heat Q is passed from one end of the plate 
eo che other. 
2. Single Plate 
Swift (1989) illustrated the basic principles of the 
thermoacoustic engine by examining a simple example of a 


single plate in a gas-filled tube such as we described 


above. In this example the acoustic and thermodynamic ef- 
fects are nearly distinct. In the absence of the plastic 
plate, the plane standing wave sustained in the tube has 
well-known acoustic properties. It's pressure and x-direc- 


tion velocity (along the length of the tube) are ziven as 
(usinz peak values): 


MEE (IN, 


and u, = i(pa/rma) costx/A) (II-2) 
where P4 is the pressure amplitude at the pressure antinode, 


the reduced wavelength A = A/2N = a/w, A is the wavelength, 


a is the speed of sound in the gas, w is the angular fre- 
quency, and fm is the mean density of the gas. The subscript 
1 indicates the first order time oscillation, and the sub- 
script m indicates the mean value. The total pressure can be 
given as 

DIS oe (II=3) 
A similar expression represents the total velocity. Follow- 
ing Swift, we will assume that the first order in the acous- 


tic amplitude is sufficient for all acoustic and thermo- 


dynamic variables (pressure, velocity, temperature, density, 


and entropy). 
The sound wave is adiabatic and has an oscillatory 
temperature given by: 


Ta = CT.B/?mce IA (II-4) 
where 8 = - (de/dT)se/Pm is the isobaric thermal expanc mom 


coefficient and Ce is the isobaric (constant pressure) heat 
capacity per unit mass. Note that T, and p, are in phase. 


For ideal gases, thermodynamics tells us that 


where Y - the ratio of isobaric to isochoric specific head 
(5/2 for monatomic gases, smaller, but greater than one, for 
other gases). Combining Equations II-4 and II-5 gives us: 
Ia/Tm = [(%-1)/¥] Pı/Pm. (II-6) 

The introduction of the plastic plate into mmens 
standing wave modifies the original, unperturbed temperature 
oscillations. This modification is due to the heat flow 
between the gas and the plate, as described previously. The 
temperature is modified in both magnitude and phase for gas 
about a thermal penetration depth away from the plate. 
According to,Swrtt- this results in two important effects: 
1) a time-ave:iage heat flow near the surface of the plate, 
along the direction of acoustic vibration, and 2) the 
generation or absorption of real acoustic power near the 
surface of the plate. In refrigerators the effect 1S an 


absorption of acoustic power. In prime movers the acoustic 
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power is generated by the temperature gradient present in 
the plate. 

There are several assumptions we will make to 
Simplify calculations, again following Swift. We assume that 
the plate is short enough compared to a reduced wavelength 
(Ax « A) and far enough from both velocity and pressure 
nodes that p, and us can be considered uniform over the 
entire plate. We assume the gas has zero viscosity, so that 
ui does not depend on y (we had already made this assumption 
by looking at u, in the x-direction only). We assume that 
the plate has a large enough heat capacity per unit area 
that its temperature does not change appreciably at the 
acoustic frequency. We assume that the plate has a mean 
temperature gradient in the x-direction v¥Tm. Finally, we 
neglect the plate's and gas' thermal conductivity in the x- 
direction. 

Applying these assumptions, we see that the mean 
fluid temperature (Tm(x)) is the same as that of the plate. 
Swift calculated the oscillating fluid temperature using 
the general equation of heat transfer. He kept only first- 
order terms, neglected thermal conduction along x, and 
applied the boundary condition 1T,(0)=0 imposed by the plate. 
The resulting equation is: 


Baty) = (ime / mee) Pa~(9Ta/m)Ua lh 1-expl-(14i)y/81]. (II-7) 


To interpret this equation, we look at it in the limit that 


2 


the gas is far enough from the plate (y”?5K) to make  negii- 


gible thermal contact with the plate. This gives: 


The first term in Equation  II-So 199 temperature 


oscillation due to the adiabatic compressions and expansions 


of the gas (see Equation II-4). The second term comes from 
the mean temperature gradient inthe gas. As the gas 
oscillates along x with displacement amplitude ua/w, the 


temperature at a given point in space oscillates by an 


amount  vT, ui/o even if the temperature oi a given piece of 


fluid remains constant. The actual temperature oscillations 
are a linear superposition of these two effects. 


The y dependent part of equation II-7 is complex. It 


approaches 1 for y » 8k, as given previously. It approaches 
zero for y « àk, where the plate imposes the condition 


Ta = 0. Most importantly, for y  5x, its magnitudes Scania 
of the order of 1, but it has a substantial imaginary part. 


This phase shift in the oscillating temperature of the 


standing wave at y Y rk, due to the thermal presence af the 
plate, is an important result because it leads directly to 
the time-averaged heat flow in the x-direction. It is the 


Same phase shift we emphasized in the description of the gas 
parcel motion. 
a. Heat Flow 
Swift argues that since we are neglecting  ordi- 


mary thermal conductivity in the Xx direction» the only way 


heat can be transported along x is by the hydrodynamic 


transport of entropy. carried by the oscillatory velocity 
ui: 


The line above the quantity s,u, indicates the time-average 
of the first-order entropy and velocity product. The heat 
flow is a second order quantity, signified by the subscript 
p. 

The total heat flow Qa along the plate, in the x- 
Enn cccron» s found by integrating d over the y-z plane: 

@> = FT of®qe dy (II-10) 
where ff is the perimeter of the plate in the y-z plane. 
EubDSpaEuting for T. and performing the integration gives: 

Q2 = -1/4 ITSk(TAB)paua(T-1) INED 
where Mr is the thermodynamically active area in a plane 
perpendicular to the longitudinal acoustic motion,  T&48 is 
the heat parameter of the gas, and T is defined as the ratio 
of actual temperature gradient to the critial temperature 
gradient  (vT/vTz-1€). The critical temperature gradient 


occurs when the temperature change along the plate just 


matches the temperature change due to the adiabatic 
compression of the gas, and no heat flows between the gas 
and the plate. It is the boundary between the refrigerator 


and prime mover functions of the thermoacoustic engines. 
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up the temperature gradient from cold to hot and wor 
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(acoustic power) is absorbed, as for a refrigerator UE 


— 


T = 1 there is no heat flow. When T > 1, heat flows down the 
temperature gradient from hot to cold and acoustic power is 
produced, as for a prime mover. 

Note that the total heat flow is proportional to 
the area lléx, and to Tm (=1 for ideal gases). It is also 
proportional to the product Pxauza, and so equals zero AS 
plate is at either a pressure node or a velocity node of the 
standing Wave. The maximum value of the product occurs 
haz vay between the nodes. Finally, the heat flow is 
proportional to the temperature gradient factor T-1. For 
l4 ^ vTaríis, T-1 > 0 and the heat flow is toward the press 
sure node, while for vT. «—vTz-43-. P1 0 and T = eee 
is away from the preszure node. If suitabie heat «=> 
at temperatures TH and Te are installed at the endz of the 
plate (with T4 - Te = 7T, Ax). this heat flow carries heat 
from ane exchanger to the other. 

The heat flow is small under ordinsmm 
circumstances. However, in closed resonators it is possible 
to achieve sound amplitudes many orders of magnitude higher 
than those of ordinary conversation. Since OS is 
proportional to p,=, and since in practical acoustic enc 
the entire cross section of the standing wave is filled with 
plates (spaced roughly 48 apart), very high "heat "UD 


and/or a large T. - Tz may be achieved. 


b 
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b. Work Flow 


The work flow (i.e. acoustic power) is given by 
the work per cycle times the rate at which that work occurs 
(the acoustic frequency f). From thermodynamics, the average 


acoustic power produced per unit volume is: 


IA e p - 9 eB p: Im[Til. [Uus 
The gas about a thermal penetration depth away 


from the plate "breathes," because of thermal expansion and 
contraction, with the right time phase with respect to 
oscillating pressure to do (or absorb) net work. This is 
exactly the same gas that we have seen 1s responsible for 


the heat flow. Gas elsewhere is ineffective in doing (or 


EHosserbing) werk. The density oscillations for y «e STER eis. 


y > 85k are in phase with the pressure oscillations, and 
hence do (or absorb) no net work. 


The total acoustic power Wa produced is found by 


integrating we over all space, as with the heat flow: 
Wo = 1/4 Mk (Y-1) (p13/Pma) (T-1) AX/A, (II-123) 
where all the terms have been defined previously. 
The acoustic power is proportional to the volume 
TSEÁX of fluid that is about a thermal penetration depth 
from the plate. It is proportional to pi, and so is quad- 


ratic in the acoustic amplitude (as was the heat flow) and 


vanishes at pressure nodes. Finally, Ws. is proportional to 
(T-1), the same temperature gradient factor as appeared in 
Q2 (Equation II-11). When vT&4 = Yrearae. (T-1) = 0, and there 
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are no temperature oscillations in the fluid other than 


those due to adiabatic processes, and no acoustic power is 
absorbed or generated. For yTm > ¥Tenge, [= MI 
tic power is produced near the plate. Whether this power 
increases the amplitude of the standing wave, is radiated 
away to infinity. is simply absorbed, or flows through an 
acoustic-to-electric transducer to generate electric power 
depends on details of the resonator, not on the plate itself 


or on the standing wave near the plate. For VIm < VTcra+, i- 


1<«< O and acoustic power is absorbed near the plate. For a 
tube without plates which has a diameter less than the 
wavelength, at constant temperature (vT = 0), this work flow 
is responsible for the ordinary thermal attenuation of a 
sound wave (Kinsler, et al.. 1982). 
E. Efficiency 

We combine Equations II-11 and II-13 to get the 

efficiency of the plate with no viscous or longitudinal 


conduction losses: 


N= W2/Q2 = (¥-1)/(TmB) (w Ax p1)/(Pma“u1).  (II-14) 
Since Uy, = uo sin x/A7M ando, "Wer toco sean 
n = (¥-1)/(T a8) (MEA Canes A (II-15) 
For X € A, n = (¥=1)7( tee) x7 R (II-16) 
and N = Ncarmor/T (Il-17) 


This efficiency approaches the Carnot efficiency 


as the power output and heat transfer rates approach zero. 


AI a A ima lar Mea leultation for the 
refrigerator mode Gf operation. where the relevant 
efficiency is the coefficient of performance, COP = Q>=/W=. 
We find that 

em cm Tea m (II-18) 

SU SNS Stack 

After developing the basic principles of the thermo- 
acoustic engine using a simplified example, Swift developed 
a more realistic model by including viscosity, longitudinal 
pacrmal conductivity, finite (instead of infinite) plate 
heat capacity, and many plates. He still made the assump- 
tions that the temperature spanned is much less than the 
absolute temperature, that the length of the plates is much 
shorter than the reduced wavelength, and that the spacing 
between the plates is greater than a few penetration depths. 
Using the equation of motion for the gas, the boun- 

dary condition that the velocity at the gas-plate interface 
1S Zero, the continuity equation of the gas, the heat flows 
in the gas and plates, and basic thermodynamic properties, 
Butt derived a wave equation for pressure P(x) in terms of 
dIm/dx, material properties and geometry. Swift also derived 
an equation for enthalpy flow along x in terms of pi{x), 
A), geometry, and material properties (including Prandtl 


number of the gas ¢ = ceu/k, where u is the viscosity, and + 


is the thermal conductivity of the gas). 


In the boundary layer approximation, for a Some 
(Ax € A) stack of plates, and neglecting viscosity, the heat 
flux and acoustic power are given by expressions very simi- 
lar to the simple results of the single plate example. 
Practical engines can be expected to have all available 
cross-sectional area filled with plates spaced about 48. 
apart. The only effect of longitudinal (1.e.. alone 
thermal conductivity is to add to the heat flow in the x 


direction by simple conduction. For gases this effect is 


normally negligible. 

The inclusion of viscosity adds considerable 
algebraic complication, anda little conceptual Gift lem 
as now ua is a function of y. 

Without viscosity, VlIz-«e 15 Be temperature 
gradient for which the temperature oscillation at a given x- 


location is zero. But in the presence of "Viscose a 


depends on x as well as y, so no well-defined vT.,.,. exists. 

To the lowest order in Ya (where gy is the fluid's 
Prandtl number), the hydrodynamic heat flow is increased by 
1/(1-48). This factor arises because the mean velocity «u,» 
underestimates the velocity with which entropy is convected 
along the stack. As we saw in the single plate calculations, 
the convective entropy transport occurs mostly at a distance 
5k from the plate. The velocity there is higher than the 


mean velocity by 1/(1-yV7). 
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The lowest order effect of viscosity on the acoustic 


power is the power dissipated by viscous shear in the fluid 
within the boundary layer, a well-known fluid mechanical 
result. 

For accurate calculations, the full general equation 
must be used; but for rough estimates, the results presented 
here are good approximations, and are much simpler to 
compute with. In addition, the results presented here are 
fairly easy to interpret physically. The expressions for 
arbitrary viscosity are complicated, but those to lowest 
order in viscosity are simple. It is easy to see the single 


plate expressions for heat flow and acoustic power about 4 


thermal penetration depth from the plates, modified as 
longitudinal thermal conductance adds to the heat flow, as 
viscous shear dissipates acoustic power, and as non-ideal 


plate properties modify the thermal boundary condition at 
the fluid-plate interface. 
4. Design Considerations 

One of the important considerations in the desien of 
a thermoacoustic refrigerator is the location of the stack 
of plates with respect to the pressure and velocity nodes 
and antinodes of the acoustic Standing wave in the resonant 
tube. Experiments performed by Wheatley, et al. (1986) 
Showed that the heat flow as a function of position fit a 
simple sine curve whose spacial period is half the wave- 


length of the acoustic standing wave. By noting how the sign 
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of the temperature difference varied with respect to the 
plate's location in the sound wave, they saw that heat 
always flows in the direction of the closest pressure anti- 
node. This effect is expected because a parcel of gas moving 
in the direction of a pressure antinode is  compressionally 
warmed and will transfer heat to the plate. A parcel moving 
toward a pressure node is cooled by expansion and will draw 
heat from the plate. At both the pressure nodes  (velocitv 
antinodes) and antinodes (velocity nodes) heat flow drops to 
zero. Thus the acoustic heat flow depends on both the acous- 
tic pressure and the fluid velocity. 


AS a plate or stack of plates is moved away from the 


pressure antinode, the temperature gradient developed 
Decomes smaller. At a quarter of a Wavelength, no gradient 
forms (the critical temperature gradient equals zero). This 


positioning effect is important in the design of a refriger- 
ator, because, together with the length of the plates, it 
places an upper limit on the maximum temperature drop 


possible across the stack. 


B. DRIVER REQUIREMENTS 


In order to drive the thermoacoustic refrigerator we 


needed a device, a transducer, that would convert the elec- 
trical energy of a battery pack to acoustic energy, as 
discussed in Subsection I.A.3. We wanted as high a pressure 


(or force) and volume velocity (or motion) aS possible to 


get the most acoustic power provided to the system. There 
were several types of electroacoustic transducers to choose 
from:  electro-static, electret, electrodynamic, magneto- 
strictive, electro-strictive, and piezoelectric. Only one of 
these, the electrodynamic transducer, can provide our system 
with both the large forces and velocities necessary. This 
transducer will be described in the next section. The 
remainder of this section describes the rest of the trans- 
ducers and explains why each was not appropriate for use in 
the thermoacoustic refrigerator. 

The electrostatic transducer consists of a pair of 
charged electrodes, or capacitor plates, one of which is 
held stationary while the other, the diaphragm, moves in 
response to electrical excitation. The electret is a kind of 
electrostatic transducer that has a diaphragm of polarized 
plastic, and therefore does not need an external voltage 
supply to the diaphragm (Kinsler, et al, 1932, p. 350). If 
an alternating voltage is applied across the plates, the 
diaphragm moves, thereby radiating an acoustic wave. 
Unfortunately there is very little force (or pressure) and 
displacement associated with this motion. This type of 
transducer may be useful as a microphone but is not an 
appropriate choice to drive the thermoacoustic refrigerator. 

In a magnetostrictive transducer a change in the 
magnetic polarization in the material causes an elastic 


Strain (Wilson, 1985, p. 3.). Although the forces are more 


than adequate the strain produces too little motion to be 
useful to us. 

In electrostrictive and piezoelectric transducers an 
externally applied electric field causes a change in the 
dielectric polarization of the material, which in WEM 
causes an elastic strain (Wilson, 1985, p. 3). Once again 
there is sufficient force but the displacement due to the 
strain is too small to be useful. 

Because these techniques have sufficient force, it would 
be possible to overcome the displacement limitation using 
mechanical leverage or other forms of mechanical impedance 
transformation, but this was not considered due ta the 


associated increase in device complexity. 


GC. ., ELECTRODYNAMIC THEORY 
1. Physics 

The electrodynamic driver is a moving coPL 
transducer which converts electrical energy to kinetic ener- 
zy (motion). It consists of a diaphragm attached to a cylin- 
drical coil of wire (voice coil) that is suspended in a 
fixed magnetic field (see Figure II-2). If an alternating 
current is supplied to the coil, the interaction of the 
current and the magnetic field will induce a force on the 
coil so that the diaphragm moves. The magnitude of the force 
is F = BlI., where B is the magnitude of the magnetic induc- 


tiong l is the length of the coil and I is the current. The 


Figure II-2. 





Schematic of an Electrodynamic Transducer 
Showing (A) Magnet, (B) Back Plate, 

(C) Pole Piece, (D) Acoustic Seal, 

(E) Diaphragm, (F) Voice Coil, 

(G) Surround (Suspension), and 

(H) Electrical Lead Connection. 


cy 


direction of the Torce given by the cross product of Gi 


current vector and the magnetic field (from Ampere's law E = 
Il X B), is everywhere orthogonal to the current and Che 
radial magnetic field. The motion of the diaphragm produces 
sound waves with a frequency equal to that of the alterna- 
ting current. Another important property of the electrodyna- 
mic transducer is that the motion of the voice coil induces 
an electromotive force (emf) which equals Bl times the 
velocity of the voice coil (emf = Blv). This "back emf" 
increases the electrical impedance of the driver and can be 
used to monitor the driver motion as discussed in Section 
EPI 

The driver is attached to the resonator via a 
reducer and bellows. The reducer is a cone-shaped piece of 
aluminum mechanically connecting the driver voice coil which 
produces the force to the bellows. which in turn is in 
contact with the gas in the resonator. The reducer is 
designed to be lightweight, rigid, and strong for a direct 
transfer of motion and energy from the moving coil to the 
bellows. The bellows is a lightweight and flexible gastight 
seal between the driver and the tube. An aluminum plate is 
used as an interface between the reducer and the bellows to 
keep the bellows face rigid. For the STAR driver the reducer 
and plate are incorporated into a single unit machined from 


aluminum bar "Stocks 
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Figure 11-23 shows a mechanical model of the 
electrodynamic transducer as a forced simple harmonic 
EEnator (SHO), which consists of a moving mass (m,,), 
stiffness (k), mechanical damping or resistance (Rm), and 
driving force (F). Included in the moving mass are the 
masses of the voice coil, reducer, stiffener, and bellows. 
The stiffness includes the suspension stiffness of the voice 
coil surround, the stiffness of the gas volume trapped 
behind the diaphragm, and the stiffness of the bellows. The 
equations governing the SHO can thus be used to determine 


the moving mass and stiffness of the driver. The primary 


relationship used is: 


or, the resonant angular frequency squared equals the 
stiffness divided by the moving mass. The use of this and 


other relationships will be discussed further in Chapter 
III, Driver Parameter Measurements. 
Ze Equivalent Circuit 
An electrical circuit can be drawn that is 
equivalent to an electromechanical system if the values 


given for Rm, Lm, and Cm in the circuit produce the same 


electrical effects as the electromechanical system itself. 
In the case of the electrodynamic driver, the equivalent 
circuit is shown in Figure II-4 with the values given by 


Kinsler, et al. (1982) as: 


Mechanical Acoustical 


(P. Glee 





e coil Ec reducer 


EZ 
susp bellows ges vol. 


Ppisurectrt-3. Electrodynamic Transducer--Forced Simple 
Harmonic Oscillator Model. 
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Figure II-4.— Equivalent Electrical Circuit for the 
Electrodynamic Transducer. 
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Rm = (BI): PaF ( II-209 
La = (BDS (I1-21) 
Cm = mo/(B1)2, (I1-22) 
where Bl is the transduction coefficient, Em LS WEE 
mechanical resistance, and k and m, are the stiffness and 


moving mass as described previously. 

The properties of the equivalent circuit allow us to 
evaluate the performance of the electrodynamic driver as if 
it were an entirely electrical system. The amount of power 
the driver generates w.1ll depend on the load it sees. Since 
one of the primary purposes of this thesis is to maximize 
the ratio of the acoutic power delivered to the electrical 
power provided (the electroacoustic efficiency), it ees 
necessary to provide a load to the system in order to 
evaluate driver performance. This load ie origina Di 
represented as a complex electrical impedence in the 
equivalent circuit. The actual value of the load will depend 
on the design of the resonator but can be characterized by 
an acoustic impedance which can be transformed ta an 
equivalent electrical load impedance. 

3. Losses 

In our search for the highest electroacoustic 
efficiency it is useful to understand where possible losses 
may occur and do what can be done to minimize them. Losses 
that are an intrinsic part of the electrical driver are eddy 


current losses in the magnetic structure generated by the 


PR 
ho 


current in the voice coil, Joule heating losses due to tne 
drive current dissipated by the voice coil electrical 
resistance, and mechanical losses due to the internal fric- 
tion of the surround and the viscosity of the gas in the 
"gap." These losses need to be taken into account. Very 
little can be done about the eddy current losses and the 
viscosity of the gas which is necessary to improve heat 
removal from the voice coil. The losses due to surround 
(suspension) mechanical resistance can be controlled by 


choice of the materials and removal of extraneous suspension 


material (see Section III.H). 


Bee COUPLED SYSTEM 
1. Computer Model 
Any discussion of the driver section of the 
thermoacoustic refrigerator is incomplete without some 
understanding of how the driver parameters are coupled to 
the resonator (acoustic load). The equivalent electro- 
mechanical circuit for the coupled system is the same as 


that shown in Figure II-4 with the acoustic load impedance: 


41 - -PoCSj[(1-j«/k)/(1-(«x/k)92) + (IIZ) 
Mos (EÉLJTSTN(KEGW+ SSiNnFéxL)coshlaeaL)TY -+ 


(Sin=(kL)cosh=(aL) + GOS eom aL) Al: 
Where e. is the density of the gas, c is the speed of sound 
(labeled "a” in previous discussions), s is the cross- 
sectional area of the resonator, k is the wave number 
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(= w/c), « - wo/(2Qc), Q is the quality factor of the resos 
nator, and L is the effective length of the resonator 
(Kinsler, et al., 1982). This circuit combines the equiva- 
lent circuit of the electrodynamic transducer with the  com- 
plex impedance load representing the resonator. One way to 
examine this coupling is through computer modeling. 


LT Michael P. Susalla, USN, developed a computer 


program which models the thermoacoustic refrigerator as the 

equivalent electrical circuit. The program models each of 

the mechanical and acoustical components as equivalent 

electrical components, as discussed above. Using the 

electrical equation from the combined circuit we determine 

the system performance as different parameters are varied. 
a. Input 


There are tuo sets of input parameters needed 


for the model: (1) driver parameters and (2) resonator 
parameters. The driver electrical parameters include the 
Voice coil inductance (LE) and resistance (RE), and the 


current (I) supplied to the voice coil (which is held at one 
ampere for all test cases). The driver mechanical parameters 
used in the model are the same equivalent resistance, induc- 
tance, and capacitance as those developed for the equivalent 
circuit. Therefore the parameter values needed for input are 
the transduction coefficient (BL), the mechanical resistance 


(RM), the moving mass (M), the suspension stiffness (SSUP), 
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the volume of the gas behind the diaphragm (VOL), and the 
cross-sectional area of this volume (A). 

The resonator parameters needed include the tube 
quality factor (Q), the speed of sound (C), the resonance 


frequency (FO), the density of the gas in the tube (RHO), 


the resonator cross-sectional area (S), and the equivalent 
length of the resonator (L). For this model the piston area 
is set equal to the resonator cross-section (S). For the 


tests reported here they were nearly equal but could differ 
sustantially in other designs. 
O AOUT puUTt 

Four output plots are generated by the graph 
program that accompanies the model: (1) voltage across the 
driver vs. frequency (since we put unit current into the 
driver, this graph can also be read as the input electrical 
impedance), (2) piston velocity vs. frequency, (3) acoustic 
power delivered to the load (resonator) vs. frequency, and 
(4) electroacoustic efficiency vs. frequency. Each of these 
plots should produce peaks at both the driver and the tube 
resonant frequencies. 

The model uses the electrical equation for the 
coupled equivalent circuit to calculate voltage directly 
(holding the current constant at one ampere) for the first 
muctput plot. It also directly calculates the force on the 
moving mass, F = BII. Maa nicals L Orec the piston velocity 


Mor the second plot is found by dividing the force by the 


mechanical impedence (Zmecn). Pressure is calculated as tube 
impedence times volume velocity, and acoustic power is 
pressure times volume velocity times the cosine of their 
phase difference (Equation I-8) for the third plot. Forum 
final plot the electrical power is calculated as voltage 
times current times the cosine of their phase difference 
(Equation I-7), and the model takes a ratio of the electri- 
cal to the acoustical powers for the electroacoustic eam 


ciency. Since the components of these powers are complex and 
have phase differences, the actual powers are calculated as 
one half the real part of [current times complex conjugate 
Or voltage], and one half the real part of [pressure times 
the complex conjugate of volume velocity]. Each of the Joere 
mentioned plots is the discussed output parameter versus 
frequency. 
2.  Assumed Typical Parameters 

Table II-1 lists the input parameters with sample 
values for the Altec 290-16K driver in air with a long 
straight resonant tube attached (see Chapter III). Sample 
outputs of the four plots using these input values are given 
in Figures II-5 through II-8. 

Figure  II-5 plots the driver voltage as a function 
of frequency. Since the current is held constant at one 
ampere, this plot actually gives us information on driver 
electrical input impedance vs. frequency. Notice that there 
are two peaks in this curve. The first peak occurs at a 


i 
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Symbol 


RM 


MABLE"II-=i 


Meaning 

Mechanical resistance 
Moving mass 

Spring constant 
Density of gas 

Sound speed 


Resonator X-sec. area 


hover transduction factor 
Resonator quality factor 
Resonator resonance freq. 
Back volume X-sec. area 


Back volume 

Voice coil inductance 
Voice coil DC resistance 
Driver current 


Resonator length 
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MODEL INPUT VALUES--EXAMPLE 


Units 
a 
Kg 
N/m 
kg/m»? 


m/s 


N/A 


DATA 


pa 
J 


ae A + —« 
AE ONU 


1.0X10-72 


17021880 Ss 


D TIS 


FREQUENCY (H2) 





VOLTAGE VS. FREQUENCY 
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© 
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Figure II-5. Computer Model Sample Output--Driver 
Voltage vs. Frequency for One Ampere 
Current into Driver. 
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PISTON VELOCITY VS. FREQUENCY 


= 
< 
i 
Y) 
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Figure II-6. 





144136 


Computer Model Sample Output--Piston 
Velocity vs. Frequency for One Ampere 
into Driver. 
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Figure II-7. Computer Model Sample Output--Acoustic 
Power vs. Frequency for One Ampere. 
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ELECTROACOUSTIC EFFICIENCY VS, FREQUENCY 





PREQUENCY (HZ) 


TIT ETT | 
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Figure II-8. 


Computer Model Sample Output-- 
pDucctnoacoust]cseErficiency vs. 
Frequency. 


S1 


lower frequency and has a much higher amplitude than the 
second. This higher amplitude peak occurs near the driver 
resonant frequency (fə = (i/2QÓ4)[k/m,]*^7 = 124 MHZ 
dominates the curve because it is the driver impedance tnat 
we are plotting, and the voice coil motion, and consequently 
the "back emf," is the largest at that frequency. The fact 
that there is a second peak near the tube resonance (FO = 
240 Hz) shows the effect of the coupling of the two parts of 
the system. 

In the plot of piston velocity vs. frequency (Figure 
II-6), the dominant peak again occurs near the driver 
resonance for the reasons presented above. The smaller 
coupled peak at the tube resonance is due to the change in 
the load impedance which occurs near the tube resonance. 

Figure II-7 shows the plot of the acoustic power vs. 
frequency. Now the dominant peak occurs at the tube 
resonance because acoustic power is primarily a function of 
the tube parameters. Here the secondary peak at the driver 
resonance is a consequence of the greatly increased piston 
velocity at the driver resonance putting energy into the 
smaller non-resonant acoustic load of the tube. 

In the plot of electroacoustic efficiency vs. 
frequency, Figure II-8, the primary peak at the tube 
resonance is sharp and the secondary peak at the driver 
resonance is small but visible. It is important to note that 


the electroacoustic efficiency peaks at a value of nearly 


60% even though the two individual resonances are separated 


by nearly an octave in frequency. The peak acoustic power 
for one amp as shown in Figure II-7 is 10.3 W. The losses 
(40%) are due primarily to Joule heating in the voice coil 
DC resistance (PR = 1/2 (I=Rec) = 5.3 W). This accounts for 
34%. The remaining loss can probably be attributed to 
mechanical resistance losses in the suspension. 


Appendix A contains a printout of the calculation 
program. Appendix B contains a printout of the graph 


program. 


III. DRIVER PARAMETER MEASUNEMENES 


Aa INTRODUCTION 

In Sections II.B and II.C we discussed the choice of the 
electrodynamic transducer to drive the thermoacoustic refri- 
gerator. Given this choice, we needed to develop a system 
for measuring the mechanical parameters of this type of 


driver ta be used in determining the electroacoustic effi- 


ciency. The relevant parameters (as discussed in Section 
Lind). cares (1) moving mass, (2) suspension Stiftinesce ee) 
transduction coefficient (Bl), and (4) mechanical  resis- 
tance. 


The driver we selected for the preliminary measurements 
is the Altec model 290-16K movinz-coil loudspeaker. There 
are several pieces of equipment needed to interface the 
driver with the resonator. Three of these parts--the reducer 
cone, stiffener, and dynamic bellows--become an integral 
part of the driver, and their effects on the driver parame- 
ters must be taken into account. We first measured the 
parameters of the driver alone, then added the other parts 
and determined the new parameters for the entire system. The 
driver-bellows combination was then fitted to one end of a 
resonant tube. The actual resonator used in the thermo- 
acoustic refrigerator will be a complex structure that is 


not easily modeled. For the purposes of the driver parameter 


measurements we used a resonant tube of uniform cross- 
section to allow easy comparison of experimental values to 
theory. The non-driver end of the closed tube housed several 
microphones. The electroacoustic efficiency of the system at 
resonance was determined from measurements of electric power 
delivered to the driver and the acoustic pressure at the 
microphone. The acoustic power dissipation was determined 
from the quality factor of the tube at resonance and the 
microphone sensitivity (the relationship between the pres- 


sure that the microphone senses and its voltage output). 


B. IMPEDANCE ANALYZER 

The primary instrument used in the testing and evalua- 
tion of the electrodynamic driver was the Hewlett Packard 
Model 4192A LF Impedance Analyzer. According to the manufac- 
turers Operation and Service Manual, the 4192A is a fully 
automatic, high performance test instrument designed to 
measure a wide range of impedance parameters as well as 
gain, phase, and group delay. Its two measurement display 
sections provide direct readout of the selected measurement 
parameters with 4.5 digit resolution along with the appro- 
priate units. In normal mode operation, the 4192A performs 
approximately five measurements per second. It also provides 
an average measurement mode (one measurement per second) for 
higher resolution and noise rejection. A high speed mode can 


perform ten measurements per second. 


The 4192A can provide measuring frequency, OSC voltage 
level, and DC bias voltage (impedance measurements only) 
equivalent to actual operating conditions. The sweep  capa- 
bility of the built-in frequency synthesizer and DC bias 
source permits quick, convenient, and accurate measurements. 
Measuring frequency or DC bias voltage can be automatically 
or manually swept in either direction. OSC level can be 
manually swept in either direction in 1lmV increments. Actual 
test voltage across or test signal current through the 


device under test is also measured. 


In amplitude/phase measurements, the 4192A can measure 
four transmission parameters - gain/loss, level, phase, and 
group delay. In impedance measurements, the  4192A can 
measure eleven impedance parameters - absolute value of 


impedance, absolute value of admittance, phase angle, resis- 
tance, reactance, conductance, susceptance, inductance, 
capacitance, dissipation factor, and quality factor eee 
impedence measurements have a basic accuracy of 0.1% and a 
resolution of 4.5 digits. Also, the circuitry of them 
provides direct and accurate impedance measurements of both 


grounded and floated devices. 


C. MECHANICAL PARAMETERS 
The most important parameters for a moving coil  trans- 
ducer are the moving mass, suspension stiffness, and the Bl 


(transduction) coefficient. Ihe movinz mass of the driver 


E 


includes the mass of the voice coil, the diaphragm, and the 
hydrodynamic mass of the air being pushed by the  diaphrasm. 
The masses of the reducer and stiffener are also important, 
as are the moving mass and suspension stiffness of the 
bellows. The moving mass and stiffness of the combined 
driver-reducer-bellows system is also important. The values 
for the combined system should equal the sum of the values 
of the individual parts. The stiffness, mass, and Bl coeffi- 
cient give us the information we need to construct the 
equivalent electrical circuit, as discussed in Subsection 
II.C.2. Also, we have to insure that the stiffness of the 
driver system is just large enough to produce a mechanical 
resonance at the tube resonant frequency, at which point the 
driver piston looks like a soft (or open) tube end. This may 
appear to be contradictory to having a velocity node of the 
standing wave at the piston, but it is not. It is how we get 
hizh efficiency. 

The important design parameters of the tube are its 
resonance frequency and quality factor. The length of the 
tube and the speed of sound in the contained gas determine 
the resonance frequency. 


Mero Dynamic mass Loading of Driver 


We used the driver's similarity to a simple harmonic 
oscillator to find the moving mass and suspension stiffness 
of the driver by taking advantage of the fact that the 


resonance frequency squared equals the stiffness divided by 


pA 


the mass. If we add mass to the driver, it should linear 


change the inverse of the frequency (or period) squared: 


T= = 1/£2 = 4n2(mo + m)/k. (IIIS 
where mo is the moving mass, m is the added mass, and k is 
the suspension stiffness. The instrument we used for this 
measurement was the HP 4192A Low Frequency Impedance 


Analyzer, described in section B. We added various masses of 
putty (one to nine grams) to the diaphragm and found the 
resonance frequency of the driver for each mass. The  putty 
mass was measured on the Sartorius model 2403 mass balance 
(serial number 1287). This balance has an accuracy of = 
0.00005 g. However, the putty couid be expected to gain or 
lose more mass than that due to handling. Based on mass 
measurements of the putty after use, we expect an 
experimental accuracy of + 0.005 g for the masses. A more 
significant uncertainty involving the masses arises from the 
inability to place the masses evenly on the voice coil or to 
produce a perfectly rigid coupling between the masses and 
the -VOICE Could It is not possible to assign a specific 
numerical percentage error to this uncertainty., but we must 
take it into account when evaluating the data associated 
with mass loading. A typical plot of the period squared vz. 
added mass in air is shown in Figure III-1. We used the 
least squares method to fit a straight line through the 
data. We determined the suspension stiffness (ROSE 


to be: 
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ALTEC 290-16K MASS LOADING 


PERIOD SQUARED (x10-6 S*S) 


Eigure III-1. 


y = 8.4947 + 1.4735x R = 0.99 





2 4 6 8 10 
ADDED MASS (g) 


Least Squares Fit of Period Squared vs. 
Added Mass to Determine Stiffness and 
Moving Mass of the Altec 290-16K Driver. 
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k 22.7 € 0.1 A OO ia (IIo 


and the moving mass (m2 2 intercep o PPIE 


Mo 


5. 8° 240i aoe (IIISSE 
for this diiver rne] 

The moving mass of the driver in vacuum is expected 
to be less than the value in air because in vacuum the 
driver is not pushing a hydrodynamic mass. However, we did 
not take any measurements of this driver in vacuum. 

2. Bellows Parameters 

The Servometer custom-designed dynamic bellows 
(shown in Figure III-2 with test stand) is used to provide a 
flexible gas seal between the driver and the resonator. 

The stiffness of the bellows was measured using two 
independent methods - static and dynamic. For the statii 
method we used Hooke's law: k = -F/x. The force was deter- 
mined by pulling a partial vacuum on the bellows and multi- 
plying the vacuum pressure (read on a manometer) by the 
surface area of the bellows (A = 1.1 ooa X10 E 
used a G. L. Collins Corporation model LMS-1i9V22-01 Linear 
Variable Differential Transformer (LVDT) (serial 274674) 
connected to an EG&G 128A Lock-in Analyzer to determine the 
displacement of the bellows under a partial vacuum. Figure 
III-3 shows the LVDT in its calibration stand. The calibra- 
tion curve for the LVDT is shown in Figure III-4. The least 
Squares fit slope of this curve gives us the calibration 


constant for the LVDT: 
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Ex Sure III-2. ` Photograph of Servometer Custom-Designed 


Dynamic Bellows--With and Without 
Test Stand. 
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Photozraph of Linear Variable 
Differential Transformer (LVDT) 
in its Calibration stand 
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LOCK-IN OUTPUT VOLTAGE (mV dc) 


PERIOD SQUARED (x10-6 S*S) 


LVDT CALIBRATION 
100 


y = 187.9945 - 378.2689x R = 1.00 





-100 
0.2 073 0.4 0.5 0.6 0.7 0.8 


LVDT DISPLACEMENT (INCHES) 


Figure III-4. Calibration Data for LVDT Voltage 
Output vs. Displacement. 


BELLOWS DYNAMIC MASS LOADING 
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y = 50.5616 + 8.6609x R = 1.00 
600 
500 


400 


300 





200 
20 30 40 50 60 70 80 
ADDED MASS (g) 


Figure III-5. Least Squares Fit of Period Squared 
vs. Added Mass for Dynamic Loading 
of Bellows. 
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dV/dx = -0.378 + 0.0005 VY/in = mm T (III-4) 
Table III-1 gives the static test data for the stiffness 
measurement. By averaging the values we obtained: 


kK = 4.3 fio 2 Xx 10m TM (lite 


TABLE III-1. BELLOWS STIFENESS TSE. 


MEASUREMENTERESUCTS 





LVDT Voltage Pressure Pressure Force (=P.a) Stiffness 

(mV ) (cm Hg) (xX10>P3) (N) (X107N/m) 
14.84 PE 397 á. 29 4.40 
1o E 4.13 4.69 4.47 
12.99 ZE DO 3.48 4.02 


Ihe dynamic measurement was a mass loading technique 
similar to that used for the driver. In this case, we 
measured the period as a function of added mass using an 
Endevco model 2225M2 accelerometer (serial FA80) and a 
Nicolet model 3091 digital storage ascilloscope ta record 
the free decay of the bellows after it had been disturbed by 
the blunt wooden end of a cotton swab. We then plotted 
period squared versus added mass and found the stiffness and 
moving mass from the slope and intercept of the least 
Squares fitted line (see Figure III-5) as discussed in the 


previous Section. This gave us values of: 
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K = 4.30 £ 0.05 X-40? N/m (III-6) 

and Mo = 5.84 + 0.05 g. (III-7) 

This value for moving mass includes the accelerometer cable 

and is not representative of the actual moving mass of the 

bellows. Subtracting the weight of the accelerometer gives 
us a value of approximately: 

mo =P POE (I11-8) 


for the bellows. We know the movinz mass has to be less than 


2.2 g, since that is the bellows' total mass--moving and 
static. The values obtained for bellows stiffness using the 
two different methods were within experimental error. We 
therefore felt confident using a value of kg - 4.5 ££ 0.1 X 


mee N/m in our later calculations. 
3. Combined System 

After we glued the reducer (m = 4.8 g), stiffener 
plate (m = 4.4 g) and bellows to the driver, we wanted to 
find the stiffness and moving mass of the combined driver- 
reducer-bellows system. We first tried a dynamic mass 
loading measurement (same procedure as with the driver 
alone), but found that the period squared versus added mass 
plot (see Figure III-6) gave unrealistic results. So instead 


we used the relationship: 


F = -kx = BII (S 

er k = -BlI/x = Bl (di/dV)(dV/dx) (III-10) 

and made a static measurement using the LVDT. We expected 
the additional equipment (reducer, stiffener, and bellows) 
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y = 1.8358 + 0.0724x R=0.99 
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Figure III-6. Least Squares Fit of Peried Squared 


vs. Added Mass for Combined Driver- 
Reducer-Bellows System. 


AM 290- 16K STATIC DISPLACEMENT 


LVDT OUTPUT VOLTAGE (mV) 





-200 -100 O 100 200 
CURRENT IN VOICE COIL (mA dc) 
Figure III-7. Driver Displacement vs. Voice Coil 


Current as Measured with LVDT/ 
Lock-in Analyzer. 
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to change the maving mass and suspension stiffness of the 
E cm but not the transduction coefficient (Bl). We applied 


a DC current to the driver and recorded the LVDT voltage 


output, then changed the current and recorded the new vol- 
tage. This process produced the values: AI = 283 * 0.5 mA 


and AV = 2.58 * 0.005 mV. We substitute these values and the 


value of Bl 21.5 + 0.5 N/A (experimental measurement of 


this quantity is discussed in the next section) into Equa- 
tion III-10 for a combined stiffness value of: 

m= S-> 2 0.1 X 10% N/m. (III-i11}? 
If we add the stiffness values for the driver and bellows 


eher a total stiffness of kK = 3.2 + 0.1 x 10% N/m 


A far 


the system is obtained. Thus the total experimental stiff- 
ness for the combined system differs by about 10% from the 
sum of the individual stiffnesses. Using the relationship 
MS. Lk/(2Rntf)^. and a measured resonance frequency of 


235 t 0.5 Hz, the combined moving mass is: 


nose ^at 5g. (IIT-120 
Epared to m> = 15.9 * 0.5 g for the individual » masses 
added together. The two types of mass measurement agree 


Euhin 22. 


fee TRANSDUCTION COEFFICIENT (BL) 
We determined the Bl coefficient of the driver using 
three independent methods. The first method was to remove 


the voice coil from the magnet and directly measure the 


Q 


magnetic induction using a Dowdy REL Industries model 912 
Guassmeter (serial 808) using a Hall effect probe madel 
912015 (serial 21122) (B = 2.2 + 0.05 Tesla). We Tien 
measured 1 (the length of the wire) by counting the number 
of turns of wire in the coil (under a microscope) ana 
multiplying this number (44) by the circumference of the 
coil (l1 = 44 X mw X diameter = 9.98 t 0.005 m). "Combine 
these two measurements gave us the value: 

Bl 022 05 070 MA» (III-13) 

In the second method we used the relationship: 

F = BlI = -kx or Bl = -kx/I = -k(dx/dV)(avW/dI). (IMSS 
We previously determined k for the driver alone (Equation 
III-2), and now used the LVDT and lock-in analyzer to find 
dx/di. The  LVDT calibration was discussed in the previous 
Section, with a calculated value of dV/dx = -14.9+0.05 V/m. 
We then measured the diaphragm displacement as a function of 
DC-current through the voice coil. This gave us a value of 
av/dl = 11.75 + 9.05 mV/A (see Figure III-7). Using the 
relation in Equation 11!-14 we obtain the value: 

Bl z 2131.2 m 09 19m (I Li =e 
In the third method we used the relationship: 
Ps  BLir= mE Or Blo= me 7 (lIIl=te5 
For this we measured the change in DC-current through the 
coil as a function of added mass, using the LVDT as a 


monitor to keep the displacement of the diaphragm constant. 
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Three different mass/current combinations (see Table IIA 
Ma duscede an average value: oí : 
BLAN A ASIN As (III-1i7)} 


Averaging the results of these three methods gave us a 


ho 


in 


value of Bl = 21.5 + 0.5 N/A with 2% deviation among the 


three methods. 


MABLE DII- 2. STATICMMEASUREMENT OF ALTEC 290-16K 


TRANSDUCTION COEFFICIENT (BL) 


ám_(8) ÁI (mA) Bl - mg/I (N/A 
132 5 64.4 ZEE 
7S 74.0 Tra 
Zov Jem 2 Dr 


E. MICROPHONE CALIBRATION 

MWeropnone calibration 1s an important part of the 
process leading to the measurement of the electroacoustic 
efficiency. A microphone is used to measure the acoustic 
pressure inside the resonator by converting the acoustic 
pressure to an output voltage. The ratio of the open circuit 
voltage to the pressure is called the open eurcugt 
microphone sensitivity. 

We used three microphones for the calibration: em —a 


General Radio (GR) model 1962-9611 half-inch electret 
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condenser microphone (serial 11811) with GR preamplifier 
model 1560-P42 (serial 4061); (2) an Endevco model S REI 
piezoresistive bridge microphone (serial 73MJ); and (3) an 
1/8 inch B&K model 4138 microphone (serial 19683). We first 
calibrated the GR using a GR type 1562-A sound-level  cali- 
brator (serial 19683, 114dB SPL, 250 Hz), and obtained) ame. 
sensitivity of Ma = 7.42 mV/Pa. This compares within about 
1% of the calibration certificate provided by the manufac- 
turer of 7.50 mV/Pa. Additional checks were made with a B&K 
pistonphone calibrator which substantiated the GR calibra- 
tion results. 

To calibrate the other twa microphones we attached the 
resonator tube to the 290-16K driver. We used an HP 33255 
Synthesizer/Function Generator to excite the driver at 240 
Ez the resonant frequency of the tube. We recorded the 
output voltages for all three microphones and used the 


comparison relationship: 


Ma/Va = M>/V> (111 wey 
to determine the sensitivities of the B&K and  Endevco. The 
B&K sensitivity was determined to be Mes = 0.462 mV/Pa. 


comparing within 2% of the manufacturer's specification 
value of 0.470 mV/Pa (adjusted to include the capacitance of 
the B&K model VAO160 1/78" to 1/4" adapter with an expected 4 
dB "line loss"). The sensitivity of the Endevco microphone 
was calculated to be 8.17 X 10-= V/Pa.”™ comparing womi 


of the manufacturer's specification of 8.09 X 10-6 V/Pa. "EG 
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all three microphones we measured sensitivity values within 


2% of the manufacturers' specifications. 


EM ELECTROACOUSTIC EFFICIENCY 

Huc cricounoaemusticwNefficleneve is the ratio of the 
acoustic power dissipated in the resonator to the electric 
power delivered to the driver. The electric power was 
determined using the relationship: 

Dae NUN Cos 8 (III-19) 
where V is the voltage across the driver terminals, I is the 
Semerent flowing in the driver, and 8 is the phase angle 
between the voltage and the current. The values for voltage, 
current and phase angle delivered to the driver were 
obtained using the HP 4192A. 

Calculation of the acoustic power was a little more 
complex. We started with the basic relationship for the 
quality factor of the tube: 

Q = 2xEs+/E.c, (TII-20) 
where Esr is the energy stored in the tube and Ere is the 
energy lost per cycle. The acoustic power generated in the 
tube is the energy lost per cycle times the number of cycles 
pen second, or frequency: 

Pac = Evcf - 2nRfEs-T/Q. (rum ou 
Following Burmaster (1985) the energy stored in a plane 
Standing sound wave is equal to the energy density 


integrated over the volume of gas in the tube (VOL): 


pa 


Est = ( PERRE ZA NOA (III-Z22g 
where P... is the root mean square acoustic pressure in the 


tube, s is the ratio of the specific heats of the gas in the 


tube (1.4 for air), and Pa is the ambient "pressure WENN 
atmosphere oor 1.013 X 105 Pa for these measurements) NN 
In ca acoustic pressure is the r.m.s. microphone "oUCtpe 
voltage divided its sensitivity (V/M) as discussed in the 


previous Section. The expression for the  acou 


tA 


tic power 
becomes: 
Pac = (N£V=VOL)/(QM=¥P.). (III-ZS 
The only remaining quantity in Equation III-23 that 
needs to be calculated separately from the efficiency 
measurement is the quality factor of the tube (Q). To 
determine Q we removed the 290-16K driver from the resanator 
and replaced it with an electret transducer to drive the 
tube. This was necessary because there were losses due to 
the 290-16K that made the tube Q appear to be lower than it 
actually is (we discovered this problem while investgating 
the reasons for the calculated efficiency to be greater than 
one). We used the B&K microphone for the pressure (converted 
to voltage) readings at the other end of the tube. We plot- 
ted voltage as a function of frequency, and determined the 
resonant frequency (f,5) and frequencies of the half power 
points (f., and fi). Using the relationship: 
Q = f./(fu-f1) (11100) 


we obtained a value of Q = 60 Bal Pop thes 


To determine the electroacoustic efficiency we returned 


the tube to the 290-16K driver (see Figure III-8) and drove 
the tube at resonance. We recorded the voltage, current, and 
phase angle for the driver and the B&K output voltage. The 
electroacoustic efficiency was then calculated using: 

Nea = Pac/Peiec. (III-25) 
Table III-3 lists the electric power, acoustic power, and 


efficiency for five test runs at driver currents ranging 


from 3-38 mA (with corresponding voltages 0.15-1.43 V). The 
electroacoustic efficiency ranged from 84.2 for the lowest 


Mieiver current to 75.9 for the highest. 


TABLBEEEDNMMLSSMLIBCSPOD.TONR ELECTROACOUSTIC 


EFFICIENCY MEASUREMENTS 


Perec (mW) Pac (mW) Mea (7%) 
OR SEO 0.47 34.2 
3 E ZZO TD 2 
1 ORG 7.79 77.9 
17m E p 
54 12 41.2 75.9 


oe COMPARISON TO COMPUTER MODEL 
The final task in the evaluation procedure for the 
electrodynamic driver is to compare the measured electro- 


acoustic efficiency to that predicted by the computer model 


Z3 





Figure III-8. Photograph of Altec 290-16K Driver 
Assembly with Resonator Tube and 
Endevco Microphone. 
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EnSscussed in Section II.D. We have already measured or 
calculated all the input values for the model except for the 


driver's mechanical resistance (R,) | 


1. Mechanical Resistance 
We can determine Rm by using the relationship: 

Ren eae Hee May, Oe: GDLI-26) 
where £f, is the resonant frequency, mj is the moving mass 
and Q is the quality factor of the driver. There are several 
methods available for determining the driver's quality 
factor. One is to plot the magnitude of the modified driver 
impedance vs. frequency and pick the half power points off 
the graph. Then, as with the tube, Q is given by: 

A La). (II1-27) 
where fo and f1 are the frequencies associated with the 
upper and lower half power points (located at [Zl]max/yZ. 
where [Z] is the magnitude of the driver impedance which is 
modified by subtracting the direct current (DC) resistance 
from the measured complex impadance). 

Another method is to use the relationship [Hofler, 
1986]: 
0 = 5$ £f, de/df. (III-29) 


where % is the phase angle between driver resistance (R) and 


reactance (X), the two components of modified impedance (Z = 
R + ix, where R is the measured resistance minus the DC 
(direct current) resistance). This modified impedance is 


Encebportional to the voice coil velocity. Ihe quantity do/df 
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phase angle versus frequency curve 


il 


is the slopcasc ecu 
evaluated at the resonant frequency (where the curve is a 
straight line). We used the HP4192A to obtain the data for 
these measurements. 

Figure III-9 shows a plot of modified impedance vs. 
frequency for the Altec 290-16K in air. From this graph we 
find f2=235.2 Hz, f.=249.6 Hz and f1=222.0 Hz. The @ for 
these values is 8.52. Figure III-10 shows a plot of phase 
angle vs. frequency for the same data set. Over this large 
range of frequencies the phase angle curve approximates an 
arctangent. At resonance, however, Where de/df is evaluate 
the curve is linear (see Figure III-11). Equation III-28 is 
the slope of the arctangent at resonance (f3). The slope of 
the plotted data is 71.8X10-* sec (dé was converted eii 
degrees to radians). The Q calculated from this method is 


8.44. These values for Q agree within 1%. From these values 


we get R42. 770g E £ 


We tried yet another method for determining driver 
Q. This third method measures the decay constant. Sia 
digital storage oscilloscope to record the Fringe down NS 
the driver. This technique takes less time to run the exper- 


iment and reduce the data using the relationship: 


Q = wt/2 or 9% =9207a Man (III-29) 
where w is the resonant angular frequency in radians/sec, T 
is the period in seconds, and v is the free decay time 
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Figure III-9. Altec 290-16K Modified Impedance 
Magnitude vs. Frequency. 


77 
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Figure III-11. Least Squares Fit of Altec 290-16K 
Modified Impedance Phase Angle as a 
Function of Frequency Near Resonance. 
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constant in seconds. This method will be discussed in 
greater detail in Chapter IV. 
2. Model Results 

Table III-4 lists the computer model input values 
for the Altec 290-16K. Figures III-12 through III-15 show 
the output plots for these values as discussed in Section 
EID. In this case the resonance peaks are much closer 
together, and the magnitude of the secondary peak is closer 
in value to the primary peak. Once again the driver reso- 
mance is at a lower frequency than the resonator resonance. 
It is interesting to note that the modeled efficiency (about 
600%) is 15-20% less than the measured efficiency reported in 
the previous Section. We believe that this results from the 
fact that for the computer model, a high tube Q means that 
Peemadriver tosses, represented by Rm. play a more important 
role in the efficiency calculation, and the results are very 
Sameer ve to small errors in R,. The primary ways for the 
driver energy to be dissipated are: (1) Joule heating due to 
DC resistance in voice coil and eddy currents induced in the 
magnet structure due to the time-varying magnetic fields 
generated by the AC voice coil current (small but not negli- 
gible): (2) mechanical resistance of driver (Rm): and (3) 
energy transferred to acoustic load. Another possible reason 
for the different efficiencies is that the efficiency is not 
independent of current hence the real physical system has 


some non-linearity, as Table III-3 clearly shows. The 


pi 


TABLE III-4. 


Symbol 


BM 


M 


Eo 
RHO 


O 


Meanins 

Mechanical resistance 
Movingz mass 

Spring constant 
Density of gas 

Sound speed 

Resonator X-sec. area 


Driver trans. factor 


Resonator quality factor 


Resonator resonance freq. 


Back volume K-sec. area 
Back volume 

Voice coil inductance 
Voice coil DC resistance 
Driver current 


Resonator length 


30 


COMPUTER MODEL INPUT=-ALTEe 


Units 


kg/s 
kg 


N/m 


K2/me 


290-T16K DRIVES 


Altec 
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mm 
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Figure III-12. 
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Altec 290-16K--Driver Voltage 


VS. Frequency for One Amp 
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PISTON VELOCITY VS. FREQUENCY 
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Figure III-13. 


Computer Model Output for 
Altec 290-16K--Piston Velocity 
VS. Frequency tor One wane: 
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Figure III-14. Computer Model Output for 
Altec 290-16K--Acoustic Power 
Vs. Frequency for One Amp. 
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-— ELECTROACOUSTIC EFFICIENCY VS. FREQUENCY 


FREQUENCY (HZ) 





308040 
0338230 
479432 
0410629 
360834 
200020 
238218 


IE 


YAMOd SINLOTIA /HAMOS OLLSNOOV 


Figure III-15. Computer Model Output for 
Altec 290-16K--Electroacoustic 
Efficiency vs. Frequency. 
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computer model is based on a linear theory and will have an 


REupdicjency which is independent of drive current. 


H. DRIVER RESONANCE ALTERATION 

The refrigerator must operate at the resonant frequency 
of the resonator to obtain the highest electroacoustic effi- 
ciency and to maintain stack position with respect to the 
standing acoustic wave (discussed in Sub-section II.A.4). It 
appears from the computer model results and our experiments 
so far that the operating frequency should also be near the 
driver's resonance to get the most power out of the driver. 
The resonant frequency of the resonator is determined by the 
effective length of the resonator and the speed of sound in 
the contained gas. The effective length for a tube of uni- 
form cross-section and closed end is the actual measured 
leneth. For the STAR resonator, however, the effective 
length is the length that would produce the measured  reso- 
nant frequency for the given sound speed if the resonator 
were a uniform cross-section closed-ended tube (see Figure 
I-2 for approximate shape of resonator). The sound speed can 
be altered by changing the percentage of xenon in the 
helium-xenon mixture, Since the sound speed depends on the 
Gensity of the mixture, and the density depends on the 
percentages of the individual gases present (see Equations 
IV-24 and IV-26). The gas mixture and effective resonator 


length will be chosen to maximize the coefficient of 


00 
un 


performance of the refrigerator. Therefore if we want to 
match the resonances, we will have to alter the resonance of 
the driver once the resonant frequency of the resonator has 
been determined. 

The resonant frequency of the driver is proportional to 
the square root of the stiffness divided by the moving mass. 
If we want to raise the resonant frequency of the driver we 
can either increase the stiffness or decrease the mass. 
Conversely, to lower the resonant frequency we can either 
decrease the stiffness or increase the mass. Generally we 
prefer to decrease stiffness and mass to get a larger 
velocity. Reduction of stiffness usually also reduces the 
mechanical resistance. 

For the Altec 290-16K attached to the uniform cross- 
section resonant tube in these preliminary experiments, the 
frequencies were already well matched (note the closeness of 


the peaks in Figures III-12 through III-14) and changing the 


driver frequency was not required. However, we know that in 
future applications we may need to alter, most probably 
lower, the driver resonant frequency. We developed an 


experimental procedure to do this. 

Adding moving mass to the driver would appear to be the 
easiest way to lower the driver resonance. Unfortunately the 
added mass would also reduce the velocity by increasing the 
mecahnical impedance above resonance. We therefore chose to 


experiment with decreasing the driver stiffness. Recall that 
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the driver-reducer-bellows system stiffness includes the 
stifffness of the bellows (4.5 X 10? N/m) and the suspension 
stiffness of the driver phenolic surround (29 X 10^? N/m). 
The driver suspension was the obvious choice far 
experimental alteration for two reasons: ates over Six 
times as large as the bellous stiffness, and (2) it would be 
very difficult for us to change the bellows stiffness 
Without affecting its performance. 

Since the 290-16K voice coil surround is made of a 
Phenolic material, it was relatively easy to cut through it 
with an ordinary scalpel. We made a series of successive 
cuts in the surround material and ran measurements of the 
driver resonant frequency after each cut. The results of 
this experiment are given in Table III-5. We plotted the 
phase angle between the driver's electromotive force (emf) 
and the voice coil displacement as a function of frequency. 
For the emf we measured the voltage across the driver termi- 
nals in channel A of the HP 4192A Impedance Analyzer and 
assumed that the emf was in phase with the driver current 
and hence force. We used the photonic sensor to measure the 
displacement by connecting its voltage output to channel B. 
The 4192A then displayed the phase angle between the two 
voltages representing force and displacement. The driver is 
at resonance when this phase angle equals 90 degrees. The 
quality factor (Q) is equal to one-half the resonant fre- 


quency times the slope of the phase angle vs. frequency 
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curve evaluated at resonance (see Equation  III-28). Using 
the previously determined values for moving mass (16 g) and 


bellows stiffness, we calculated the suspension stiffness: 


ae = Kae Der. Kei Lows (TII-30) 
where Keoe = 4n fao nn We calculated the mechanical 
resistance  (R,) using equation III-26. We stopped cüt EES 


the surround when its stiffness was approximately equal to 
that of the bellows and therefore the total stiffness was no 


longer suspension dominated. 


TABLE -III-S mANRESULTS TOR 


DRIVER RESONANT FREQUENCY ALTERATION 


Percent Surround Resonant Quality Suspension Mechanical 


Remaining Frequency Factor Stiffness Resistance 
(approx. (Hz) (410N/m) kg/s) 
100 ZUM) 9455 S2 2.889 
1007 SIX SLICES) po RS 8.09 28.8 2.83 
ao Spo Se by, 24a aep s 
50 122.5 8405 16.3 2148 
20 119.5 GOES 4.52 o 


Despite the anomoly at the 80% data point, we feel that 
this experiment shows that partial removal of the surround 
can reduce stiffness, resonant frequency, and mechanical 


resistance. It appears from the above data that the flexing 
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ef the phenclic surround is responsible for the majority of 


the mechanical losses in the driver. 
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IV. STAR DRIVER OBARADMTBPBHESS 


À. INTRODUCTION 

This chapter discusses the selection of the electro- 
dynamic driver, describes the parameter measurements on the 
chosen driver, and estimates the electroacoustic efficiency 
of the STAR using the computer model described in Section 
Dn 


ls Interface uith NASA GAS Cannister 


In addition to the maximum weight of 200 pounds 
discussed in Section I.B, restrictions due to the Get  Awav 
Special (GAS) cannister dimensions zive us a maximum payload 
height (28.25 inches) and diameter (19.75 inches). In our 
alloweu height we must fit two battery layers, the driver 
with its housing, the resonator assembly, and its associated 
vacuum space and thermal insulation. The battery layers are 
each about four inches high, for a total battery height of 
eight inches. This leaves about 20 inches for the driver and 
resonator assemblies (see Figure I-4 for approximate 
configuration). The electronics package and gas reservoir 
should fit in the area radially surrounding the resonator, 
and will therefore not add to the height. Although the 
diameter of the cannister volume is 19.75 inches, the driver 
housing must attach to the standard GAS can lid at a 9.5 


inch diameter bolt circle, otherwise we would need a lid 
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specially designed for our project. As far as weight is 
concerned, approximately every pound we save elsewhere 
allows us to carry an extra battery cell, which results in 
an extra ten watt-hours of available energy. 
2. Maximizing Electroacoustic Efficiency 
In addition to the space and weight restrictions, we 


want to maximize the electroacoustic efficiency of the 


driver. We can do this by obtaining the maximum value for 
mie combination of parameters  Bl/R,m5. A larger Bl will 
give us a larger driving force (F=Bli), which in turn gives 
us a larger acoustic volume velocity, one of the components 
of acoustic power. For a given value of magnetic induction 
(B), a longer voice coil (achieved by having a larger voice 
coil diameter with the same number of turns. or the same 
diameter and more turns) will result in a larger Bl. This 
longer 1 would also result in a greater voice coil 
electrical resistance, so the two effects cancel. We also 


want to minimize the moving mass (m2) to allow greater 
velocities and minimize the driver's mechanical resistance 


CA) in order to minimize mechanical losses in the driver. 


3. Driver Selection 
We had two off-the-shelf drivers available for con- 
sideration. The Altec 260-16K was discussed in Chapter III. 
The other driver uas a JBL 2445J. Both drivers were designed 
for applications somewhat different from the STAR. Figure 


IV-1 shows a side view of the 2445J driver. The hatched 





Figure IV-1. Schematic Drawing of JBL 2445J Driver. 
Hatched Area Shows Unnecessary Throat 
Section. 
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wea. known as the throat of the driver, is considered 
unnecessary for our application, and adds a lot of extra 
weizht to the driver. The Altec 290-16K has a similar excess 
throat volume. We contacted JBL and asked if they would 
custom-design a driver which would eliminate this excess 
weight. Fortunately for our project, JBL was in the process 
of designing a new driver, the 24503, which uses neodymium- 
iron-boron magnets instead of the traditional ceramic. They 
agreed to custom-design a version of the 24503 for our use. 


The resulting driver (hereafter referred to as NIB) is shown 


in Figure IV-2. The NIB driver accepts the same voice coils 
(four inch diameter, aluminum wire, and titanium suspen- 
sion/diaphragm) as the 2445J (ceramic magnet) driver. SEL 


claims that it should also have the same properties, 
Mme luiding Bl. 

In the remainder of this chapter we describe the 
parameter tests on the NIB driver. We then compare the NIB 
Emyer to the 290=16K driver... The results of the 24453 
parameter measurements are given for comparison. Finaliy we 
present results generated by the computer model for the NIB 


driver given various acoustic load impedances. 


B. NIB DRIVER PARAMETERS 


For the measurement of the NIB driver parameters, W 


a 


used the techniques we developed in Chapter III for th 


(D 


measurement of the 290-16K driver parameters. 
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1. Moving Mass and Suspension Stiffness 


The moving mass and suspension stiffness were 
calculated using the mass loading technique. We placed small 
masses (0.5 to 3.5 grams) of black putty on top of the 
voice coil in half-gram increments. For each added mass the 
resonant frequency of the driver was found using the HP 
4192A LP Impedance Analyzer (see Section III.B for a 
description of this instrument). Figure IV-3 shows a sample 
plot of period squared vs. added mass for the NIB driver in 
air. The stiffness was found by calculating four pi squared 
divided by the slope. The moving mass is equal to the 
intercept divided by the slope. For the NIB in air these 


values are: 


= eee. 1X 10 ke /s= (IV-1) 
and mg = WB £ 0.1 E, ve) 
with a resonant frequency of 564 Hz. We repeated the 


measurement in a vacuum chamber and found: 
K= 8.6 + 0.1 X 10% kg/s2 (IV-3) 
and Mg O e. (IV-4) 
with a resonant frequency of 833 Hz (see Figure IV-4). 

The moving mass is greater in air because of the 
driver having to move the air itself. The discrepancy in the 
stiffness is not as easy to explain. To determine the actual 
value of k, we used the static mass loading method. This 


gave us a value of: 


NIB MASS LOADING - AIR 


y = 3.1368 + 0.4022x R=1.00 
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Figure IV-3. Least Squares Fit of Period Squared vs. 
Added Mass for NIB Driver in Air. 
NIB MASS LOADING - VACUUM 
2.6 | | P 
y = 1.4403 + 0.46x R=0.99 
gp 2.4 
e 
e 
A 22 
[C 
S 
3 2.0 
0 
e 
O 1.8 
CE 
Lu 
CE 
1.6 
1.4 
O 1 2 3 
ADDED MASS (G) 
Figure IV-4. Least Squares Fit of Period Squared vs. 


Added Mass for NIB Driver in Vacuum. 
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E OA ALO? ke/s=. (IV-S) 

(The static loading measurements will be discussed in tne 
next Sub-section. ) 

This last value for suspension stiffness differs by 

about 40% from the first two. We did not have sufficient 


time to test this driver until we were sure we had the 


"correct" value for the stiffness. We recommend that the 
follow-on students do sa. Lacking the correct value, we ran 
the computer model for two sets of data: one set represen- 


ting the "high" k (9.0 X 10% kg/s?) and one using the "low" 
k (6.4 X 10% kg/s2). 
EXC ransdauction Coefficient (Bl) 
The transduction coefficient (Bl) was calculated 
using two methods. The first method was to measure B (magne- 
tic induction in the voice coil gap) and 1 (length of the 


voice coil) directly. We measured B with the Dowty RFL 


Industries model 912 Guassmeter (serial 308) using a Hall 
effect probe model 912015 (serial 21122). The maximum 
reading using this probe was 18.7 kGauss = 1.87 Tesla (T). 
This value differs by less than one percent from the value 
stamped on the driver by the manufacturer: 1.858 T. We will 
therfore use a value of B=1.86 + 0.01 T for this driver. We 
tren counted the turns of aluminum wire on the voice coil (n 
= 30 | turns) and calculated 1 using the relationship: 


LO A (IV-6) 
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For the 9.98 * 0.005 cm diameter coil we calculated a value 
of 1 = 9.4 730.3 m. Multiplying these two values together 
gives us for the transduction coefficient: 

Bl  1779S^* 9637» (IV-7) 


The other method used to determine Bl used the 


relationship: 

F = Bll = -kx or BL = =K P (IV-8) 
where I is the direct current applied to the coil, k is the 
spring constant measured earlier, and x is the displacement 


of the voice coil producea by the current. For this measure- 


ment on the 290-16K driver we had used a Linear Variable 


Differential Transformer (LVDT). For the NIB driver, 
however, we used Mechanical Technology Incorporated's MTI- 
1000 photonic sensor (serial 38010818) with pluz-in 


cartridge model MTI-3812 (serial 1366) and probe model 125R. 
Figure  IV-5 shows the photonic sensor with the NIB driver. 
One major advantage of the photonic sensor over the LVDT is 
that the photonic does not need to make physical contact 
With the surface it's measuring (it uses light) whereas the 
LVDT does. With the Altec 290-16K's rigid phenolic dia- 
phragm, the contact with the LVDT probably did not affect 
the measurement. The NIB diaphragm, however, is made of thin 
titanium and is quite flexible. It would have been  impos- 


3ible to get an accurate readinz of the displacement of this 


g9 


diaphragm using the LVDT. We calibrated the photonic sensor 


using calibration fixture KD-CH-III A. The output voltage 
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Figure IV-5. MTI-1000 Photonic Sensor with NIB Driver. 
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from the sensor was read on a HP 3478A multimeter (serial 
A110136). We then plotted DC voltage vs. displacement irn 
millimeters. The resulting graph is shown in Figure  IV-eo. 
The slope of this line gives us the calibration constant: 

dV/dx = -2.086 +4 0.001 nME (IV-9) 
This value is within two percent of the calibration value 
provided by the manufacturer of dV/dx = 2.045 V/mmn See 
slope). 

The next step was to run DC current through the coil 
and record the displacement of the diaphragm (via a voltage 
reading from the photonic sensor). For the voltage supply we 
used the HP 467A power amplifier (serial 444-00243). We 
found the current by measuring the voltage (with the HP 
3476A digital multimeter, serial 292162) across a ten ohm 
(102) resistor placed in series with the circuit, and 
dividing the voltage by ten ohms. The resulting graph is 
shown in Figure IV-7. The slope of the line gives us: 

dlzdV 2 2.10-*-0. 009 2742 (IV-10) 

Combining the values in Equations IV-9 and IV-10 with the 

previously calculated high end value for the suspension 
stiffness k - 9.2 5 0-6 Xx 107 Ec cM E e. 

Bl = k(dx/dV) (dV/dil) = 2125 $T (IV-11) 

This 15% difference between the values of Bl using 
the two independent methods is what prompted us to test for 
Bl and k using the static loading method. Our primary 


concern is that we would expect the direct measurements of B 
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PHOTONIC CALIBRATION 
y = 11.0854 - 2.0353x R= 1.00 


DC VOLTAGE (V) 
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Figure IV-6. Calibration Data for Photonic Sensor 
Voltage Output vs. Displacement. 
[ NIB - PHOTONIC CURRENT CALIBRATION 
= - 15.2167 + 2.1004x R = 1.00 
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6.9 720 TA 7272 ves 7.4 755 


VOLTAGE (V) 


Figure IV-7. Least Squares Fit of DC Current Through 
Voice Coil vs. Photonic Voltage Output. 


TON 


and 1 to give us the highest possible value for Bl. This is 


because the maximum magnetic field may not be acting on the 
full length of the coil. We did not expect the value of Bl 
calculated using the spring constant to be 15% higher than 
the direct measurement value. 

We were able to obtain static loading values for 
both Bl and k uSing the same equipment setup and the 
Formulas 

Bl = Amg/AI and k = - Amg/Ax (IV-12) 
For these measurements the voice coil is driven with a DC 


current. We then add previously measured masses (putty and 


ul 


oldering wire) to the voice coil, and measure the resulting 


displacement of the diaphragm using the photonic sensor. To 


Q 


measure Bl we determine how much current is needed to bring 
the voice coil back to its null displacement (original 
output voltage reading from the photonic sensor) after the 
mass is added. We measure k by keeping the current constant 
and measuring the difference in photonic output voltage 
(multiplied by the calibration constant to get displacement) 
resulting from the added mass. Table IV-1 lists the results 
of these measurements. 

The bracketed value ineach k and Bl pair was 
calculated using the measured value in the pair and the 
relationship: 


Bil =) ee (IV 


Lo» 


From these measurements we get the value k = 6.4 * 0.2 kg/sz 
as reported earlier, and the transduction coefficient: 

Bl = 13580 * 0.2 N/A. (IV-14) 
For the model runs we used a value of Bl = 19.0 for the hizh 


k and Bl = 15.0 for the low k. 


DCER. STATIC LOADING OF NIB DRIVER TO DETERMINE 


SUSeeNSION STIPPNESS AND TRANSDUCTION COEFFICIENT 


åm (8) AI (mA) AV (V) k (X10* ke/s?) Bl (N/A) 
i. 1 a — 5.2 6.7 [15.3] 
n. 1 -- 5.5 6.4 [14.6] 
34.7 ze -- [6.4] 14.6 
34.7 22.3 -- [6.7] cu 
34.7 -- .288 6.1 -- 


(front slope) 


3. Mechanical Resistance 


The next parameter we measured was the mechanical 


resistance (R.) using the relationship: 


Rm = 2nfono/O. (IV-15) 


We attempted to find the driver quality factor (Q) using the 
Same modified impedance technique described in Section 


III.G. As can be seen in Figure IV-8, the graph of the 


modified impedance magnitude vs. frequency is not a smooth 
resonance curve, as it was for the 290-16K (see Figure III- 
9). Bn fact, it looks like the superposition of two 


IZ| (OHMS) 


NIB IMPEDANCE MAGNITUDE 





480 990 


Figure IV-8. 


520 540 560 580 
FREQUENCY (HZ) 


NIB Modified Impedance Magnitude 
VS. Frequency. 
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different impedance curves, the first one shorter and fatter 
and the second one taller and thinner. It was impossible to 
accurately measure half-power points from this curve, So we 
chose the phase angle method. (It is interesting to note 
that the JBL 2445J (ceramic) driver did have a smooth modi- 


fied impedance curve similar to that of the 290-16K.) Figure 


IV-9 shows the graph of phase angle vs. frequency over a 
wide range of frequencies. Although this curve generally 
approximates an arctangent, the segment of the line immedi- 


ately adjacent to the resonant frequency (assumed to be the 
frequency where @=0(, in this case 534.5 t 0,5 Hz) is linear 
(for a frequency range of about ten Hz). This segment is 
EEtted in Figure IV-10. The slope of the line is de/df = - 
42.0 * 0.5 X 10-?s (after converting degrees to radians). As 
fmm section I1.G, Q is given by: 
CA Mer dE) = 11.2 2 0.2. (IV-16) 
Mis sives a value for Ra (using the moving mass in air 
Mo = 7.8 + 0.1 g): 
Ree 2234 270805 ke/s. (IV-17) 
We decided that a more accurate measurement of Q and 
Ba. could be found by making the Q measurement in vacuum. For 
this we used the open circuit free decay method. This method 
consists of sending a signal to the driver at resonance with 
an HP 3314A function generator (serial A342044), then 
simultaneously switching off the signal and triggering a 


Nicolet model 310 storage oscilloscope (serial 88004641) -tO 
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Figure IV-9. NIB Modified Impedance Phase Angle 
as a Function of Frequency--Full 
Frequency Range. 
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Figure IV-10. NIB Modified Impedance Phase Angle 
as a Function of Frequency--Linear 
Segment Near Resonance. 
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record the free decay. A sample of the output from this 
process is shown in Figure IV-11. The equation for a simple 
harmonic oscillator experiencing free decay is: 

A(t) = Ao exp (-t/t) sinl(w# t+), (IV-18) 
where A(t) is the amplitude as a function of time, Aj is the 
initial (maximum) amplitude, t is the time (= nT, or the 


number of periods times the period), t+ is the free decay 


constant (= QT/rx, see Equation III-29), wo is the angular 
resonant frequency, and # is the phase angle. We recorded 
the peak amplitude as a function of period number, and 


plotted the natural logarithm of the peak amplitude vs. 
period number. The equation representing this line is: 
A A = MT IT. (IV-19) 
The slope of this line is t/t = nT/t = nrx/Q. Therefore Q can 
be found by dividing n by the slope of the line. Figure IV- 


12 shows a semi-log plot for the NIB driver in air. The 


value of Q for this plot is 21.5, which is almost twice the 
value we found using the modified impedance phase angle 
measurements (Q - 11.2). Obviously the air strongly affects 
the measurements. The scattering of the data points in this 
figure is probably due to the same phenomena that causes the 
double peaks in the impedance curve. Fizure IV-13 shows the 
semi-log plot of the amplitude vs. period number for vacuum. 
This curve is much smoother and, we assume, more reliable. 


Tne Q for this data is: 


Q = x/slope = 19.8 + 0.05. (IV-20) 
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Sample of Free Decay Output from 
Nicolet 310 Storage Oscilloscope. 


Figure IV-11. 
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NIB FREE DECAY - AIR 
y = 3.5792 -0.1464x R = 1.00 
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Semi-log Plot of Amplitude 
vs. Period Number for NIB 
Free Decay in Air. 


Figure IV-12. 
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y = 4.5606 - 0.1585x R= 1.00 
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Figure IV-13. Semi-log Plot of Amplitude 
vs. Period Number for NIB 


Free Decay in Vacuum. 
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We can also find the period from the free decay curve (and 
thus the resonant frequency as the inverse of the period). 


We use the oscilloscope to measure the time for ten cycles 


and divide this number by ten. For the vacuum run the 
resonant frequency was 839 t 0.5 Hz. These values (and the 
moving mass - 3.2 g in vacuum) give us the mechanical resis- 


bance in vacuum: 
Rm = 2NfoM9/Q = 0.83 2 0.01 ke]. (IV-21) 
This is the value we used for the model. 
Table IV-2 gives a comparison of applicable 
parameters for the three drivers. Notice that the parameter 


Bl/Rmácrmo is almost three times greater for chemia m TE 


oy 


than the 290-16K. This is the parameter that we said (Sub- 
section IV.A.2 we wanted to maximize to get the highest 
thermoacoustic efficiency. Thus the driver that zives us the 
best height and weizht characteristics (see photograph of 
the three drivers in Fizure IV-14) also gives us the best 


electroacoustic efficiency. 


Cs MODELED EFFICIENCY 

In order to model the electroacoustic efficiency of the 
driver, we need to know what acoustic load (impedance) it 
will see. The final resonator and stack composition will not 
have been determined by the time this thesis is completed. 


We therefore assumed a range of possible values for the 
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TABLE IV-2. COMPARISON OF PARAMETERS FOR 
THE THREE DRIVERS 

Property (units) Altec 290 JBL 2445] SIN TB 
Moving mass. m, (kg) 0.0151 0.015 DES 
mering constant (N/m) 3.52x10* "OST pero 
Magnetic Induction, 

B (Tesla) 2.1 1.97 IS 
Voice Coil 

Diameter (cm) PAD 9.98 9.98 
Bl (N/A) Dco 1375 TANS 
Voice coil DC 

resistance (a 10.6 Pal 7.0 
Mechanical 

resistance (kg/s) AS, “a. 0.82 
merece coil 

inductance (H) B O0xiOT- LO ie IO 

(at 20 kHz) 
Mass (kg) i uS OE 
Height (in) OT) 52.5 2 
Magnetic enerzy 

density (J/cm?) eos 0 ae 
EN Rina Mo LO 1230 1410 

(m/Coul-kg) 
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acoustic load based on the experimental data acquired to 
date. 

The model parameters that relate to the acoustic 
load are the density of the gas (RHO), the speed of sound in 
the gas (C), the resonator quality factor (Q), the resonator 
resonant frequency (FO), and the resonator length (L). The 
model was designed for a closed-ended tude of uniform cross- 
section. The STAR resonator has an open tube (with a bulb 
attached) and nonuniform cross-section. The problem of 
applying the existing computer model to the STAR resonator 
was solved by scaling the STAR parameters to make it look 
like a closed-ended uniform cross-section tube. 

1. Input Values for Pure Helium 

For pure helium, the experimental data available was 
the range of resonator quality factor values (10-20), and 
the speed of sound and gas density at room temperature 
(about 300K) and ten atmospheres (Creer = 1024 m/s. 
fer = 1.618 kg/m). Also available were four resonant fre- 
quencies, two for a short resonator at warm (580 Hz) and 
cold (520 Hz) temperatures, and two for a long resonator at 
Warm (520 Hz) and cold (470 Hz) temperatures. Using this 
information, we first calculated the effective lengths of 
the resonators at room temperature as if they were  closed- 
ended, uniform cross-section tubes. For this we used the 


relationship: 


AMAS COL. Lis» = c/2f, (IV-22) 


is 


where A is the wavelength, L is the length of tne tube NM 
the speed of sound, and f is the resonant frequency. This 
calculation gave us effective lengths for the resonators: 
Le = 0.883m and Li = 0.985m. For the computer runs cR LLLI 
convenient value Of Larr = 0.90m. We then calculated the 
speed of sound for the operating temperature of approxi- 
mately 250K using the relationship: 
C(T) = “Grr I teal F (IV-23 

This gives us a speed of sound Coso = 925 m/s. If we use 
this sound speed and the effective resonator lengths to 
estimate the resonant frequencies for 250K, we get fy, = 529 
Hz and fi = 474 Hz, which agree with the measured values. We 
calculated the gas density at the operating temperature 
using the relationship: 

ATA DA a mm (IV-24) 
where y is the ratio of the specific heats (1.667 for 
helium) and P is the pressure. 

Table IV-3 lists the input parameters for the pure 
helium computer runs. Note that for the computer model, the 
moving mass must take the mass of the reducer cone and 
bellows into account (see Figure IV-15). The total moving 
mass is expected to be approximately 15 z. 

2. Input Values for Helium-Xenon 

There was somewhat less data available for the 

helium-xenon mixture (12.5% Xe). We had the sound speed at 


300K and ten atmospheres (Crr = 472 m/s). We also hau 
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NAME E OMBUTER MODEL INPUT-- 
DEBEMPESDUEBGUITUSHELIOMCATITIEN ATMOSPHERES 
AND TEMPERATURE = 250K 


RS 


Symbol Meaning umts NIB 
RM Mechanical resistance kg/s NUES 
M Moving mass kg DOTES 
SSUP Spring constant N/m 6.6, 9.0 X10=- 
RHO Density of gas kg/m? 1.94 
C Sound speed m/s 935.0 
= Resonator X-sec. area m= A A >? 
BL Driver transduction factor N/A SS SO 
Q Resonator quality factor TOO Sado 
FO Resonator resonance freq. Hz 519.4 
A Bocla volumne x Sec. arca mz 4.42X10-9 
VOL Back volume me 2 ZO MO € 
LE voice coil inductance H 1c 15 DIS 
RE Voice coil DC resistance ái 7.6 
J Driver cürrent A TO 
L Resonator length m 0-90 
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Figure IV-15. Photograph of JBC NIBET eE 
with (a) Reducer and (b) Bellows. 
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resonant frequencies for the same four conditions as for the 


Sere helium: short-warm (26S Hz), shert-cold (230 Hz), long- 


warm (235 Hz), and long-cold (215 Hz). We used Equation 
Eoo to find the effective resonator lengths (La = 0O.891m 
and L, = 1.004m). However, once again we used a convenient 
effective length of 0.90 m for the computer runs. We used 


Equation IV-23 to find the sound speed for He-xe at the 


operating temperature, Caso = 431 m/s. The estimated reso- 
nant frequencies for the cold temperature were f.-242 Hz and 
f,=215 Hz, once again comparing favorably with the measured 
Values. To find the density of the mixture at the operating 
temperature we used a form of the ideal gas law: 

EM 7 (IV-25) 
where M is the atomic mass of the mixture, R is the 
universal gas constant, and P and T are the pressure and 
temperature. We calculated the atomic mass of the mixture 
using the equation: 

ieee = Mia t Mae. CIV- 6) 


where the M's are atomic masses (4 for helium and 131.3 for 
xenon) and the x's are percentages of the gases in the 
mixtures. This equation gives us an atomic mass for the 
12.5% Xe mixture of 19.9 g/mol, and a density of 9.70 kg/m>œ 
at 250K and ten atmospheres. Table IV-4 lists the computer 
input values for He-Xe, including a range of Q's from 5 to 
ED. (We are not sure what to expect for the actual value of 


Q. since a neu stack will also be used.) 


EA 


Symbol 


RM 
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SUP 


RHO 


AQ 


ft 


BE 


FO 


VOL 


LE 


RE 


IV-4. COMPUTER MODEL ANC 


AND TEMPERATURE = 


Meaning 
Mechanical resistance 
Moving mass 
Spring constant 
Density of gas 
Sound speed 

¿or X-sec. 


Reso- area 


Driver trans. factor 


Resonator qualitv factor 


Resonator resonance freq. 


Back volume X-sec. area 
Back volume 

Voice coil inductance 
Voice coil DC resistance 
Driver current 


Resonator length 
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3. Computer Model Results 

Figures IV-16 through IV-19 are sample output plots 
for pure helium with Bl = 15 N/A, k = 6.4 X105 pe/s2,. and 
Q = 20. Notice that the peak values for the driver and 
resonator are almost matched (the driver resonance is at the 
lower frequency). The thermoacoustic efficiency for this set 
of values is about 70.4%. Figures IV-20 through IV-23 show 
the corresponding plots for the helium-xenon mixture. This 
time the peaks are not matched (the driver peak is off the 
scale) and the electroacoustic efficiency is about 16.6%. 
The primary reason for the difference between the two gases 
is the lower sound speed, and therefore lower resonant fre- 
quency, in the He-Xe mixture. 

Table IV-5 lists the maximum value of 
electroacoustic efficiency for each of the test runs. The 
consistently low efficiencies for the He-Xe emphasize the 
importance of matching the two resonant peaks. We recommend 
that the driver frequency be adjusted to match the resonator 
operating frequency. This can best be done by reducing the 


stiffness of the titanium surround (see Section III.H). 
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VOLTAGE VS. FREQUENCY 
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Figure IV-16. Computer Model Output for NIB with 
Helium (Bl = 15, Q = 20) -Driver 
Voltage vs. Frequency for One Amp. 
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PISTON VELOCITY VS. FREQUENCY 
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Figure IV-17. Computer Model Output for NIB with 
Helium (Bl = 15, @ = 20)--Piston 
Velocity vs. Frequency for One Amp. 
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Figure IV-18. Computer Model Output for NIB uith 
Helium (Bl = 15, Q = 20)--ACoustic 
Power vs. Frequency for One Amp. 
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Figure IV-19. Computer Model Output for NIB with Helium 
(Bl = 15, @ = 20)--Electroacoustic 
Efficiency vs. Frequency. 
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Figure IV-20. 
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Computer Model Output for NIB with Helium- 
Xenon (Bl = 15, Q = 20)--Driver Voltage 
vs. Frequency for One Amp. 
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PISTON VELOCITY VS. FREQUENCY 
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Figure IV-21. Computer Model Output for NIB with Helium- 
Xenon (Bl = 15, Q = 20)--Piston Velocity 
vs. Frequency for One Amp. 
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Figure IV-22. Computer Model Output for NIB with Heliums 
xenon (Bl = 15, Q2 20)-SACOUStl1CWMESU D 
vs. Frequency for One Amp. 
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Figure IV-23. Computer Model Output for NIB with Helium- 
xenon (BI = T5. Q - 20)--Electroacoustic 
Efficiency vs. Frequency. 
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TABLE IV-5 

He (Bl=15) He (Bl=19) 
le LOM 
71.6 74.3 
704 76.2 
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Vac eer VEER CONSTRUCTION 


A. INTRODUCTION 
This chapter describes the construction of the STAR 
driver. A large part of this thesis involved the design of 
the parts necessary to integrate the NIB driver with the 
resonator and GAS cannister, as illustrated in Figure I-4. 
The major design work was for the driver housing. The 
requirements for the driver housing are: 

1 - lightweight, 

2 - strong (pressurized to 210 psi), 

3 - interface to GAS canister lid (12 10-32 bolt holes 
centered at 9.5 inches diameter) , 

4 - interface with resonator (preferably interchangeable 
with earlier test resonators), 

5 - helium leak-tight feed-through ports for electrical 
leads (driver electric power, accelerometer and 
microphone outputs), and 

6 - hold NIB magnet structure securely. 

Figure V-1i shows a cross-sectional drawing of the driver 
housing (not including any of the feed-through ports, which 
will be discussed later). Distinction between parts is shown 
by different hatching. Part A is the NIB magnet structure 
which was discussed in Chapter IV. Part B is the main body 


of the driver housing, machined out of 6061-T6 aluminum. 
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Figure V-1. Cross-sectional Schematic of the STAR Driver 
Housing Showing (A) NIB Driver, (B) Driver 
Housing Main Body, (C) Aluminum Pusher Plate, 
(D) Pressure Lid, (E) Leak Test Lid, 
(F) Resonator Vacuum Housing, and (G) Bolt 
Holes for Interface with Resonator. 
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Part C is a 6061-T6 aluminum pusher plate used to hold the 
magnet structure securely in the body using six 10/32 bolts. 
Part D is the pressure lid. sealed a 1/16th inch diameter 
NE c-ring and held in piace with six 10/32 bolts. Part E 
represents part of the leak test lid (see Figure C-1). When 
completed this will be replaced by the GAS can lid provided 
by NASA. The driver housing must match up with the CAS can 
Eu 9.5 inch diameter 10-32 bolt circle. Part F is the 
resonator section housing. The arrows labeled G point to the 


bolt holes which form the driver-resonator interface. 


m DESIGN CONSIDERATIONS 

The driver-resonator system will be pressurized to ten 
atmospheres with a helium-xenon mixture (12.5% xenon). 
NASA's regulations say that it must be pressure tested to 
1.5 times the operating pressure (NASA, 1986a), or 15 atmo- 
mares (z1.52X10e Pa = 220 psi). An added complication of 
working with helium (or helium mixtures) is that it diffuses 
through almost all materials except metals. As a result, we 
are not able to use rubber o-rings to seal metal parts 
together, but must use lead o-rings instead. 

1. Electrical Feed-Throuzhs 

NASA's regulations prohibit us from leaking gases 

into space. primarily to avoid contamination of other exper- 
iments and for the safety of the Shuttle and astronauts. 


Tits regulation is important not only for the major o-ring 
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seals, but also for the microphone and accelerometer outputs 
to the frequency control electronics (see Figure I-5) and 
connecting the audio amplifier to the driver. There are 
three individual wires and three co-axial connectors that 
must penetrate the driver housing without allowing leakage. 
The other housing penetrations are a DC pressure gage for 
the driver housing and a fill line from a helium-xenon 
bottle which will stabilize the pressure in the system 
against temperature changes in the environment and those 
generated by the refrigerator. 

Based on prior experience in other applications, we 
decided that Emerson & Cuming's STYCAST*T"^ 2850FT epoxy 
(Canton, MA 02021) would be the best choice for providing 
the electrically insulating seal for the wires in the 


housing. STYCAST«rm» 2850FT is a quartz-filled epoxy formu- 


lation with remarkably good over-all general properties. La 


addition to having excellent electrical grade insulation 
properties, good adhesion, and excellent resistance to chem- 
icals and solvents, 28S0FT has unusually high thermal con- 
ductivity and low thermal expansion. This combination of 
properties makes this epoxy valuable for applications where 
electrical insulation and mechanical protection must be 
maintained while coping with heat transfer considerations. 
One definition for a sealed container is a helium 
leak rate of less than 10- standard cc's per second (NASA, 


1986b). To test the epoxy we designed and built a test 
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fixture to be used with the Veeco MS-9 helium leak detector. 
A cross-sectional diagram of the test fixture is shown in 
Figure V-2. Part A is a copper tube with 1-1/8" outer diam- 
eter to interface with the leak detector.  Soldered to the 
copper tube is a piece of solid brass stock about one inch 
thick (B) with a 0.052 inch diameter hole bored through the 
center. Epoxied into the hole is a 0.025 inch diameter bare 
copper wire (C). This design is motivated by the desire toa 
provide as long a travel path and as small a cross-section 
as possible for the helium gas molecules. Part Dis a hollow 
aluminum housing designed to hold pressurized helium. We 
decided that a rubber o-ring should hold the helium suffi- 
ciently well for a short term leak test. Part Eis a 1/8" 
NPT (National Pipe Thread) hole for connecting to a helium 
fill line. A photograph of the test fixture with and without 
the housing is shown in Figure V-3. 

The copper tube and brass stock were soldered toze- 
ther. We boiled the assembly to remove excess  particulates 
and rinsed it with alcohol. We then etched the hole in the 
brass and the copper wire with nitric acid diluted with 
distilled water (about 20% acid) to provide a clean, rough 
Surface for the epoxy. Before etching, cellophane tape was 
placed on the portion of the brass stock where the o-ring 
would be sealing (after making sure that the tape would not 
dissolve in the acid). We then rinsed the assembly with 


distilled water and dried it, mixed the epoxy, and inserted 


Figure V-2. 





Cross-sectional Schematic of the Helium 
Leak Detector Test Fixture Showing 

(A) Copper Tube for Leak Detector Interface, 
(B) Solid Brass Stock, (C) Bare Copper Wire, 
(D) Aluminum Housing, and (E) Connection for 
Helium Fill Line. 





Figure V-3. Photograph of the Helium Leak Detector 
Test Fixture. 
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the wire into the hole in the brass stock, making sure the 
epoxy filled the gap. 

The following procedure was used to run the leak 
detector tests: first zero the meter, then do a test run on 
the sensitivity calibrator. The sensitivity calibrator used 
was a Veeco type SC-4 serial number 3179 with a leak rate of 
©.5X107-% std cc/sec (standard cubic centimeters per second). 


We would then remove the calibrator, recheck the meter for 


N 


era, attach the test fixture, read the meter, remove the 
fixture, and recheck for zero. 


The test fixture was initially pressurized to 115 


psi with helium. We tried to test it immediately but found 


ou 


tnat the detector would not draw encugh of a vacuum to check 


for helium. We thought that there might be outgassing from 
flux or something in the tube. 3o we pulled a vacuum on the 
tube overnight and the next day it worked fine. We tested 


the fixture on the first, third, and sixth days after MEM 
Deen pressurized. The numerical results are listed in Tatle 
V-1. Even after we had waited six days for the helium to 
diffuse through the epoxy, the fixture leakage rate waS an 
order of magnitude less than the calibrator. The leakage 
rate was therefore four orders of magnitude less than that 
allowed by NASA to be considered leak-tight (1xX10-4 «std 
cc/sec). After the third leak test we checked the pressure 
in the aluminum chamber and found it ta be within a ten 


percent of the original pressure. We therefore concluded 


that the bare wire-epoxy method with a length/diameter ratio 


of about 20 could be used for the driver housing penetrating 


electrical feed-throughs. 


TABLE V=1. HE TUM] LEAK DETECTOR TEST RESULTS 


GONECODDEROWERESEBOXY FIXTURE. 


Number days after null meter  calibrator copper wire feed 
pressurization reading reading -through reading 
x Q nur 0202 
Es O eZ QUOS 
5 0 tas l ora 


Figure V-4 shows a close-up of the feed-through plug 
for the wire leads. It will house a total of five wires 
(with space available for a sixth if necessary). Two of the 
leads will be for the voltage supply to the voice coil. 
These leads will consist of bare copper wires epoxied into 
the small holes seen in the figure. Each wire will then be 
soldered to a 6/232 phillip-head screw with a 1-72 blind- 
tapped hole in it. The screw will be glued into an 8/22 


hole, forming a mechanically stable, demountable, electrical 


terminal. The accelerometer and microphone leads will also 
have bare copper wire epoxied into the small hole, and will 
Ee connected to Microdot''"* coaxial fittings (Malco 





Figure V-4. Driver Housing Feed-Through Plug 
for the Electrical Leads. 


Pibe OS 0001001, south Passadena, CA 91030). 
Figure V-S shows an example of each of these leads. 
2. Measurement Devices 

The microphone and accelerometer used by STAR are 
amiar to those used by Hofler (1986 and 1988). En facet. 
the resonator end of the driver housing was designed to 
accept Hofler's resonator, calibration volume, and instru- 
mentation for testing. The microphone used to measure the 
acoustic pressure in the resonator is a small quartz-crystal 
(Valpey-Fisher, Hopkinton, MA 01748) or PZT dynamic pressure 
transducer (Piezoelectric Products lot number 4-9-84, 
FO13184, Metuchen, NJ). It is mounted in the driver housing 
where the housing connects to the resonator flange. The 
microphone will be calibrated by comparison to a standard 
transducer, which will have been previously calibrated via a 
reciprocity technique. The microphone will be connected to a 
miniature FET follower circuit. The FET buffer (Eletec model 
304 impedance converter, Daytona Beach, FL 32020) converts 
the microphone's high impedance output signal to a low 
impedance to reduce sensitivity to electrical noise and 
maintain the output voltage and phase while driving the 
cables necessary to reach the signal processing electronics. 
The buffer's DC power supply and output cables are complete- 
ly  shielded through the driver housing to the buffer (see 
Figures  C-3 and C-7 for location of microphone and FET in 


driver housinz). 


PS 


(a) 





(D) 
Figure V-5. Example of Wire Lead Connectors 


for Feed-Through Plug: (a) Single 
Wire and (b) Co-axial Cable. 
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The Endevco model 22 Picomin E Eme no meters 


attached to the side of the bellows stiffener facing the 
driver. The integrated accelerometer signal will be used as 
a measure of the bellows velocity, which will be calibrated 
by driving a small cavity of calculable acoustic impedance. 
foes the product of the outputs from the microphone and 
integrated accelerometer can be combined to yield the acous- 
tic power and their relative phase can be used to keep the 
driver frequency at the tube resonance. The operation of the 


refrigerator at resonance is essential bath to optimize th 


(D 


ma E 
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acoustic power delivered to the "stack" and to  guarant 
that the "stack'" is always located at the proper position in 
the standing wave as discussed in Sub-section I1.A.4. 

The DC pressure transducer is an  Onmega-*T"? Px-20 


designed for high pressure applications requirin 


m 
Fa 
O) 
J 

in 
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stability. Figure V-6 shows a drawing of the pressure trans- 
ducer and its driver housing fitting. 

The helium port in the housing is provided ta 
empty and fill the driver and resonator and so that a stable 
pressure may be maintained in the driver and resonator 
through the ballast volume provided by the reservoir (see 
Figure 1-4). 

There is a 1/32 inch diameter hole filled with a ten 
thousandths inch inner diameter half inch long capillary in 


the driver housing where it connects to the resonator. The 


Meet lion ot the capillary is to equilibriate the static 
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DRIVER 


TRANSDUCER HOUS I NG 
FITTING 





Figure V-6. DC Pressure Transducer and Driver 
Housing Plus Fitting. 


pressure between the driver and resonator during filling 
non-operating times while maintaining a pressure seal 
acoustic operating frequencies. 

A complete set of dimensioned drawings for 
driver housing and its auxiliary equipment is given 
Appendix C. A list of materials for the driver portion 


the STAR is given in Appendix D. 
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VI. CONCLUSION AND RECOMMENDA TIHE 


The objective of the STAR project is to test and space 
qualify a new continuous cycle cryogenic  refrigsensbsm 
system for cooling of sensors and electronics which is based 
upon the newly-discovered thermoacoustic heat pumping 


effect. The new refrigerator has no sliding seals, a cycle 


frequency of about 300 Hz, and uses acoustic resonance to 
enhance overall efficiency. This thesis is concerned speci- 
fically with the design and testing of the electrodynamic 
transducer which is responsible for the electro-acoustic 
power conversion. A computer model of the driver/resonator 
system was presented along with the techniques for measure- 
ment of the electrical and mechanical parameters used as 
input for the model. 

A final driver design (including dimensional drawings) 
Vel Ua Zine a modified  JBL 2490 neodymium-iron-boron 
compression driver and associated leak-tight electrical 
feed-throughs, microphone, accelerometer, pressure gage, 
pressure housing, and resonator interface is provided. 

Since we only received the NIB driver in late December 
1987, there wasn't time to perform enough parameter 
measurements to obtain accurate values for the suspension 


stiffness (k) and transduction coefficien = NE The next 
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EM -:d-urt-esheurd find these values for use in the 
model and efficiency calculations. 

The resonator experiments (not reported here) indicate 
that  Helium-Xenon may be the preferred gas mixture for the 
STAR. If this is the case, the NIB driver resonant frequency 
will have to be lowered for it to match the operating 
(resonant) frequency of the resonator. Dno could be 
accomplished by adding moving mass to the voice coil- 
reducer-bellows assembly, but that would decrease the 
velocity and lower the electroacoustic efficiency. The 
preferred option for driver resonance adjustment is tc 
decrease the suspension stiffness by cutting away part of 
the surround (see Section III.H). Since the NIB voice coii 
has a titanium surround and diaphragm this operation will 
have to be performed carefully. The titanium diaphragm 
(dome) will also have to be removed. (The dome's function is 
performed by the reducer and bellows.) 

The task for the next group of students is to assemble 
the NIB driver with its housing and accessories (feed- 
throughs, microphone, accelerometer, etc.) and perform 
efficiency measurements on the ground similar to those that 
Will be performed in space. This series of measurements will 
closely resemble those performed on the Altec 290-16K driver 


as described in Chapter III. 
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APPENDIX A. COMPUTER MODET- CATCUERATTON ESSE 





PROGRAM FRIDGE 
INTRODUCTION: 


THE PURPOSE OF THIS PROGRAM WRITTEN BY LT.  MICHARISN 
SUSALLA IS TO CALCULATE THE DATA POINTS AND PLOT SRE 
PERFORMANCE CHARACTERISTICS OF A CONCEPTUAL MODEL OF A 
THERMOACOUSTIC HEAT ENGINE. THE COMMENTS INCORPORATED Se 
THIS PROGRAM AND THE COMMENTS INCORPORATED IN THE PLOT 
PROGRAM "GRAFRI FORTRAN A1", WILL ALLOW THE USER TO  PRO- 
DUCE FOUR GRAPHICAL PLOTS WHICH WILL CLOSELY APPROAT 
THE PERFORMANCE OF AN ACTUAL THERMOACOUSTIC HEAT ENGINE. 

THE GOAL OF THE PROGRAM/SIMULATIO" "S TO ALLOWS 
OPERATOR TO VARY DIFFERENT PARAMETERS TO ACHIEVE WHAT IS 
BELIEVED TO BE AN OPTIMUM CONDITION WHERE THE ENCED 
DRIVER RESONANCE IS MATCHED WITH THE TUBE RESONANCE. 


GENERAL INSTRUCTIONS: 


THIS PROGRAM IS SET UP TO PROVIDE FOUR PLOTS: VOR NM 
VERSUS FREQUENCY, PISTON VELOCITY VERSUS FREQUENCY e 
TIC POWER VERSUS FREQUENCY, AND ELECTROACOUSTIC ERE k E as 
VERSUS FREQUENCY. IN GENERAL THE PROGRAM HAS HARD TINE 
DATA SUCH THAT ALL YOU MUST DO TO GET A SET OF PLOTS 
COMPILE "GRAFRI FORTRAN Ai" IN FORTVS AND THEN RUN "FRIDGE 
FORTRAN Al" IN WF77 FROM FLIST CN THE NPS MAINFRAME USING 
A TEKTRONICS 618 TERMINAL. THIS ACTION WILL GIVE YOU AR 
OF PLOTS FOR THE INPUT DATA. TO TRY OTHER DATA YOU MUSAS 
TO XKEDIT ON THE "FRIDGE FORTRAN Ai” PROGRAM AND VARY THE 
INPUT DATA. AS LONG AS YOUR CHANGES ARE RESTRICTED TOMUS 
"FRIDGE" PROGRAM THE ONLY ACTION NECESSARY TO PRODUGES TE 
PLOTS IS TO RUN THE PROGRAM IN WF77. IF YOU FIND TITANES 
SARY TO CHANGE THE PLOTTING PROGRAM,  "GRAFRI" . FOR WEA 
EVER REASON, IT MUST BE RECOMPILED IN FORTVS PRIOR CIO 
NEAT ATTEMPT TO PLOT OR THE CHANGES WILL NOT BE SEEMS 
THE MACHINE. 


Ep SEC IE 

A. DEFINITION OF VARIABEES- 

ALL OF THE VARIABLES HAVE BEEN ESTABLISHED AS REA 

OR COMPLEX QUANTITIES BASED ON THE REQUIREMENTS TO 
ARITHMETIC  ALGORITHOM OF THE CONCEPTUAL NOE. EACH OF 
THE VARIABLES WILL NOW BE DEFENEOS 

LCOIL-THE IMAGINARY PART OF THE IMPEDANCE OF THE VON 
COIL OF THE DRIVER, AND IS CALCCULCATEDET OE PE 

K=THE PROPAGATION CONSTANT OF THE WAVES IN THE TUBE 
AND TS GIVEN MB ne. 


pur 


(Ut 


de 


Ca 
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IMZM=THE IMAGINARY PART OF THE MECHANICAL IMPEDANCE OF 
THE TUBE/DRIVER COMBINATION. IT IS GIVEN BY A RATHER LONG 
DIIDNEESSMPSEOUEEEGUATITONCIN-LHE CALCULATION SECTION OF THE 
PROGRAM. 

MAGU=THE MAGNITUDE OF THE COMPLEX QUANTITY “U”, IN 
METERS ORC. 

IMZT=THE IMAGINARY PART OF THE TUBE MECHANICAL IMPE- 
Paes aNDetls GIVEN BY SUBTRACTING THE DRIVER PORTIONS FROM 
TMAM. 

MET PERMOV INC MASS TOE THE SYSTEM AND IS THE SUM OF THE 
VOICE COIL AND THE DIAPHRAGM MASSES OF THE DRIVER, IN 
KILOGRAMS. 

E CDLBRESEDCRISORSNEIB.TUBE IN METERS. 

Ee hoe nibyueC TANCE SOF ITHE VOTCE COPE, IN HENRYS. 

oN Cae CURRENT  FROMEIHE SOURGE, IN AMPS. 

MAGV=THE MAGNITUDE OF THE COMPLEX QUANTITY "v", IN 
VOLTS. 

ZM=THE MECHANICAL IMPEDANCE OF THE TUBE/DRIVER SYSTEM, 
eNDPIT IS MADE UP OF THE REAL PART, REZM, AND THE IMAGI- 
WAR PAR IMZM. 

ZE=THE ELECTRICAL EQUIVALENT OF THE MECHANICAL IMPE- 
Pech Oh Teese TUBE/DRIVER SYSTEM. IT IS GIVEN BY TRANSFOR- 
ENDE MESACRMOSSC.ASMTRANSPORMER USING THE BL COEFFICIENT OF 
MEE DRIVER. 

OCE E COMECE ECECTRICAL [IMPEDANCE OF THE VOICE 
COTL. 

ZIOT=THE TOTAL ELECTRICAL IMPEDANCE OF THE SYSTEM, 
COMPLEX . 

v= pie COMPEEA VOLTAGE OF THE SOURCE NECESSSARY TO 
PRODUCE THE SELECTED INPUT CURRENT. 

UNUM BSOCIBIMOP THE PISTON, (DIAPHRAGM). OF THE 
DREIVERFIN M/S. 

ZTUBE=THE COMPLEX MECHANICAL IMPEDANCE OF THE TUBE 
IESELP. 

EENHMCOMEIESBOLPRESSURE OF THE SYSTEM. IN NEWTONS/M**2. 

PI=A CONSTANT FOR CONVERSIONS TO AND FROM ANGULAR 
UNTTS. 

RM=THE MECHANICAL RESISTANCE OF THE MOVING MASS OF THE 
SS TEM. 

Deo AN TCAL SIIFPENESS OF THE SYSTEM DUE TO THE 
DRIVER SUSPENSION. 

RHO=THE DENSITY OF THE GAS IN THE TUBE AT A SPECIFIED 
MEAN TEMPERATURE. 

Pee oe ee OR SOUND FOR THE SPECIFIED GAS AND 
TEMPERATURE. 

S= THE CROSS SECTIONAL AREA OF THE TUBE. 

BS CUE RENTES TEN TO OR STAE DRIVER THAT IS BASED ON THE 
ENEBDY AIRMCAP FEWX DENSITY OF THE-DRIVER AND THE LENGTH 
CERTS VOICE COIL. 

Q=THE QUALITY FACTOR OF THE TUBE. 

FO=THE RESONANCE FREQUENCY CF THE TUBE. 

A=THE CROSS SECTIONAL AREA OF THE BACK VOLUME OF THE 
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ORIVERE 

VOL=THE VOLUME OF THE BACK VOCUMETOTTTAEE 2. 

RG=THE INTERNAL ELECTRICAL RESISTANCE OP TRAEN 
SOURCE 

RE=THE ELECTRICAL RESISTANCE OF THE VOlGl we eine 

SVOL=THE MECHANICAL STIFFNESS OF THE sYSTEMSDUE ae 
BACK VOLUME "OF THE DRI Vane 

WO=THE RESONANCE FREQUENCY OF THE TUBE IN ANGULAR 
UNITS,  RAD/SEE.. 

ALPHA=A CHARACTERISTIC QUANTITY OF THE TUBE THAT 
FUNCTION OF THE RESONANT FREQUENCY, THE O AND THE RES 
SOUND. 

F=FREGUENCY, IN HERTZ. 

W=THE FREQUENCY CONVERTED TONRAD D SECA 

FORCE=THE FORCE ON THE MOVING MASS DUE TO THE Tj 
ACTION OF THE VOICE COIL MAGNETIC FIELD AND THE "DR 
MAGNET: 

REZT=THE REAL PART OF THE COMPLEX MECHANICAL 
IMPEDANCE. 

PU=THE PRESSURE TIMES THE VOLUME VELOCITY All 
POWER. 

PWR=THE ELECTRICAL POWER NECESSARY TO DRIVE THE EAS 
AL IHES PECCP EDOTNEUTEC CURRENTE 

PUPWR=A POWER RATIO, ACOUSTIC TO ELECTRICAL. AN 
CIENCIA 


B.DATA DESCRIBRTTON: 

FOR THE SELECTED INPUT DATA THAT IS “HARD INP Sea 
FOLLOWING VALUES WERE EMPIRICALLY DETERMINED: M C Or 
A, VOL, BL, 5, SSUP, AND RE. RHO AND C WERETSECEC b PE 
HELIUM GAS AT 10 ATMOSPHERES AND A MEAN TEMPERATURE OFE EEE 
DEGREES ESI. RM WAS CALCULATED FROM A THEORETICAL 
EQUATION USING RHO,C S AND THE DRIVER RESONANT FREQUENCY. 
RG WAS SET TO ZERO TO MODEL AN IDEAL POWER SOURCE. "CES 
SET TO ZERO FOR EASE OF CALCULATION, BUT WAS TESTED BURG 
RUNNING TO VERIFY THAT IT HAD LITTLE IF ANY EFFECT RE 
INCCUDED -AFA CURRENT I, OF ONE AMPERE WAS SELECTED AHORA 
EASE OF ANALYSIS. 

SINCE YOU MAY WANT TO CHANGE MY INPUT DATA I WILL PUT 
A COPY IN THE COMMENTS SO THAT YOU CAN RETURN TOS 
DESIRED. MY ANALYSIS HAS SHOWN THE BEST BACK VOLUME TCE 
9.50E-8, AND THE BEST RM TO BE 1.535 KG/S. THIS VACO E EE 
RM WAS DETERMINED AFTER REVIEW OF PLOTS SHOWED HIGHER THAN 
EXPECTED VOLTAGES AND IT WAS DETERMINED BY MEASUREMEN Gia 
THE LABORATORY. 


DATA RM.M,SSUP,RHO,C/2.5S23E-3,0.0045,2.6E4,1.68,1000/ 
DATA S,BL,Q,.FO,A/1.521E-3,8.15 7800 mr SF 
DATA VOL.RG,LE.RE,I1.,L/2.715E-6.0 0 SmE 
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ou D qusc 


Cre 


poo tee hem MAGUIMZT .M,L,LE.1,MAGV 
SO leer om oe, ceOll, ZIO@r,Vv,U.ZTUBE.P 
Pil= 22141592654 


DAA he Mo oor RHO. EEES E= oo tia, 2 ofa. 1000 / 
DATACS BC O, FO., DIN TE-S S 15. 800, EY 
DATA AAA 6,0,0,8,1 D oz 


SVOL ATAT RHOFC*C)/VOL 
WO=2*PI*FO 
ALPHA=WO/ (2*Q*C) 


MAECUEATION SECTION" 

BIENUNEUVIMESPOEINESOCREATEB ANC OUTPUT FILE FOR THE DATA 
Mmoetels NECESSARY 10 PRODUCE THE DESIRED PLOTS, AND PLACES 
HEADINGS ON THE DATA AS APPROPRIATE. 


OPEN (UNIT=69, FILE='FRIDGE OUTPUT A1',IOSTAT-IOS) 


WRITE (eo, *)~ MASS VOLUME ' 

WRITE(69,1)M.VOL 

ince (GO, * ) ” FREQUENCY ENSIS SEOSEBRMBSOSCUENSIT VEL PWR RATIC 
s VOLTAGE' 


ISBEOBMATCELS.4,E13.4) 


THE FOLLOWING DO LOOP PROVIDES FOR THE CALCULATION OF 
BENE EB HUNDRED DATA POINTS FOR EACH OF THE DESIRED PLOTS, 
USING A FREQUENCY RANGE OF ONE QUARTER TUBE RESONANCE TO 
THREE HALVES TUBE RESONANCE. 


BO 10 J=1,300 
F=(F0/4)+((J-1)*(F0/240)) 


PEE THREE CINES CALCULATE VARIOUS QUANTITIES. 
EXPLAINED ABOVE, THAT ARE FUNCTIONS OF FREQUENCY, FOR USE 
IN FUTURE CALCULATIONS. 


Wz2*PI*F 
KzW/C 
LCOIL=W*LE 


Ar ESPN TAE INES CALGULATE THE REAL AND IMAGINARY 
PARTS OF THE TUBE/DRIVER MECHANICAL IMPEDANCE AND COMBINE 
MEM INTO A COMPLEX QUANTITY, ZM. 


eevee © > ALPHA~COS(K*L)*SIN(K*L))+(K*RHO*C*S 
3 B UMDAEPHASTOSCOSHCALPHA*EL)))/OK*(1-CALPHA/K)**2) *(C(CC( 
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* SIN(K*L))**2) * ((COSH (ALPHA*L) )**2) )+(( (CCS ay ee 
SINE CAL PRA 2a A 
IMZM= (W*M)-(SSUP/W) -(SVOL/W)+( ( (RHO*C*S5*ALPHA* SINE (Ae 
E *COSH( ALPHA*L) )-(K*RHO*C*S*COS(K*L)*SIN(K* 107 Ge 
$ /K)**2)*(CCCSINCK*L)) **2) * (C CCOSHCALPHA* L)O 7E V 
$ x*2 J= (SINH(ACEHA CN A ERR 
LM=CMPEX (RE METZA 


C THE NEXT THREE LINES CONVERT ZM TO AN ELECTR 
C QUANTITY AND THEN COMBINE ALL OF THE ELECTRICAL IMPEDANGiI 
C JO GET A TOTAL IMEEDANCES 


ZEE M 
ZCOILSCMPERSORESICOTSIE) 
ZIOT=RGFZCOIL+ZE 


THE NEXT SIXTEEN LINES CALCULATE THE VALUES CLON 
PLOTTED, AND PUT THEM IN THE OUTPUT FILE, AND [THEN CSS 
THETETCES 


Pad o 


Vd mp 
FORCE=BL* 1 
U=FORGH 221 
MAGU=CABS(U) 
REZT=REZM-RM 
IMZT=IMZM-(W*M)+(SSUP/W) + (SVOL/W) 
¿LIB CM (FELT Zi 
ESAS 
PU=O 2S (REAL P* CONNIE (i 
PWR =O 7S" (REALL CONIG ys 
PUPWR=PU/PWR 
MAGV=CABS(V) 
WRITE (69,100) F ,MAGU PU, PUPWRaAey 
GO FORMAT ( E13 .5,E13.5, E13. 5 Ele) 
10 CONTINUE 
CLOSE (og) 


PIOPLPINGSESSEGTTON: 
THE NEXT TWO LINES INTERFACE WITH MY PLOTTING PROGRESS 


"GRAFRI FROTRAN A1", AND THE DISSPLA ROUTINE TO PLOT THE 
DATA. 


ooo 


PAUSE "PRESS ENTER TO START PEOTPINGOBOSBNUUM 
CALL EXCMS( EXECUODISSBULASCROHLTM 

SLOP 

END 


SO 


440404400 040000000000 (CORSO 


OQO C 


APPENDIX B. COMPUTER MODEL GRAPH PROGRAM 


PROGRAM GRAFRI 


PO NOC HAEL EP. SUSALLA'S PLOTTING PROGRAM FOR 
USE IN CONJUNCTION WITH "FRIDGE FORTRAN Ai” IN "THE THER- 
meAcOUSTIC HEAT ENGINE SIMULALTION”’. THIS PROGRAM MUST BE 
COMPILED IN FORTVS PRIOR TO RUNNING "FRIDGE FORTRAN Al” IN 
WF77.THIS PROGRAM IS FOR USE IN DISSPLA AND CALLS VARIOUS 
MISS EA” SUBROUTINES. THESE VARIOUS SUBROUTINES ARE USED 
MORE THAN ONCE IN THE PROGRAM BUT WILL ONLY BE DESCRIBED 
ONCE EACH. IN ORDER TO USE DISSPLA YOU MUST HAVE ARRAYS 
Pies DATA POINTS TO BE PLOTTED. MY ARRAYS ARE ALL DIMEN- 
ENMONEDO TO 300 POINTS SINCE THAT IS THE NUMBER CALCULATED 
BY "FRIDGE FORTRAN Al”. YOU SHOULD BE AWARE THAT IF YOU 
CHANGE THE NUMBER OF POINTS TO BE CALCULATED, AND THE 
MUMBER IS GREATER THAN 300, YOU MUST REDIMENSION THE 
ARRAYS. MOR HSRACTERSSTRINGS ARE USED IN THIS PROGRAM. 
NEW1$ IS USED IN RENAMING FILES TO THE FORMAT NECESSARY 
Pee Oloor LA, AND DUMMYS IS USED TO READ THE HEADINGS OUT 
Side DATA FILE. ALL OF THE DATA TO BE PLOTTED IS IN THE 
Meme OF REAL NUMBERS, sO ALL OF THE VARIABLES ARE MADE 
REAL. (D IS USED FOR THE FREQUENCY DATA, G IS USED FOR THE 
TON VELOCITY DATA, H IS USED FOR THE ACOUSTIC POWER 
DATA. O IS USED FOR THE ELECTROACOUSTIC EFFICIENCY DATA 
AND, T IS USED FOR THE VOLTAGE DATA. THE XMAX AND YMAX 
VARIABLES ARE INITIALIZED TO ZERO AND ARE USED TO CALCU- 
LATE THE MAX VALUE OF EACH OF THE PARAMETERS AS THEY ARE 
READ INTO THE ARRAYS. THIS COMPLETES THE INTRODUCTION AND 
SARTAIN S THE FIRST FIVE LINES OF THE PROGRAM. 


CMM ENS TONTO 3000) G C300) ,H{ 300) ,O(300) ,T(300) 
CHARACTER*80 NEW1$,DUMMYS$ 

Paes O,G,H,Ogl ,M 

DATA XMAX,YMAX1,YMAXZ,YMAX3,YMAX4/0.,0.,0.,0..,0./ 


THE NEXT SEVEN LINES OF THE PROGRAM RENAME MY DATA 
EN DIO THE FORMAT NECESSARY FOR USE IN DISSPLA, OPEN THE 
ieee READ THE MASS AND VOLUME OUT OF THE FILE AND, READ 
Pee dEADINGS OUT OF THE FILE. 


NEW1$-'RENAME '//' FRIDGE OUTPUT A1'//' FILE XYZ1 A1' 
CALL EXCMS(NEWi1$) 

SERNOUOQOSUBILEZ' XYZ1') 

READ (69,'(A80)') DUMMY$ 


mS 


CY OO) GO) 


ACI EOS 


REAL (697 1) Me ver 
READ (69, (A80 DUAL 
1 FORMAT (ELS 4 1c) 


THE FOLLOWING DO LOOP READS THE PLOTTING DATA INTs 
ABOVE SPECIFIED ARRAYS FOR USE IN THE DISSPLA SUBROUG ae 
AND ALSO CALCULATES THE MAXIMUM VALUE OF EACH OF SG 
ARRAYS FOR SETTING UP THE SCALES OF TRE Pik@ is fe THEN 
THE DATA FILES CLOSED THESE VALUES ARE LEFT AS TER EEES 
NUMBERS VICE INTEGERS BECAUSE THE DATA FROM THE TEPERT 
CAN VARY SIGNIFICANTLY WITH DIFFERENT INPUT VALUES, AND IF 
THE MAX VALUE IS LESS THAN 10 AN INTEGER REQUIREMENT GIVES 
a PLOTTING ERR. 


DO 20 J=1,300 
READ (69,15,END=21)D(J) G( IIA O 


I FORMAT(ET13.5.5813.5,E13.5 ET NN MEE) 
IP (DO TOS GO INA SUPPE 
XMAX=D(J) 

ENDIF 
IF(S(IJGEGT. MAS LI) TEEN 
YMAX1-2G(J) 

ENDIE 

IPN E UMAXS TEEN 
YMAX2=H(1) 

ENDIF 

IPC OCT GT YMAX:ITHEN 
YMAX3=0(3J) 

ENDIF 

Pett CI) eleva xa Tee 
YMAX4=T (J) 

ENDIF 


20 CONTINUE 
21) CLOSE (es) 


THE NEAT THREE LINES TAKE CARE OF ANOTHER ODDITY See 
THE DATA, IN THE CASE OF THE VOLTAGBAPIESHM DEPENDING ON 
THE INPUT DATA THE VOLTAGE WILL LOCK IN AT 8 VOS 
EXPECTED, BUT 3 THIS IS ALSO THE MAX VALUE. SOR NICE 
STRAIGHT LINE AT THE TOP OF THE PLOTOSMATICHIME THE NEL 
THREE LINES INCREASE THE MAX VALUE TO 9 VOLTS TO MAKES 
PLOTTING CASPESMORE CLEAR: 


IFCYMASA4 LE. 9 0)TDENSIN 
YMAX4=9.0 
ENDIF 


o Gn EQ qp. 


(SOON eoo TODICONXO CC eo Oc 


THE NEXT EIGHT LINES CALCULATE THE INCREMENTAL STEPS 
ROR T EACH AIS OF THE VARIOUS PLOTS. THER AIS IS THE 
AOR AEB IN INTEGER FORMAT FOR CLARITY. 
EACH AXIS IS SET UP FOR TEN INCREMENTS. Mame CAPA FILE Is 
THEN RENAMED TO THE DESIRED FORMAT FOR MY STORAGE. 


INC=INT (XMAX/10. ) 
YINCI=YMAA1/10. 
YINC2=YMAX2/10. 
YINC3=YMAX3/10. 
YINC4=YMAX4/10. 
MAX=INT(XMAX) 
NEW1$2'RENAME '//' 
CALL EXCMS(NEW1$) 


PA LL EEDSNBOE OUTPUT AI 


PES. ¿io THE MENESES ET ING SECTION. Ere bers THE 
INITIAL SET UP OF ALL OF THE CALLS THAT CAN BE COMMON AND 
SPOTS TAE VOLTAGE VERSUS PREQUENCY PLOT. «ABL TERG1S SETS 
UP THE REQUIRED MACHINE FOR THE PLOTTING AND LINKS YOU TO 


fea PLOTTER: C UEM MENMOUDEMSSSSEDEIOMAGNISY THE PLOT FOR 


Posed R PICTURE. Gane Geen UTS THE AXIS NUMBERS 
Pee I NDICULAR TO THE PLOT AXIS TO PREVENT RUNNING INTO 
EACH OTHER. SADA SE MS UP ITHE PAGET SIZE EN. INCHES. 


CALL AREA2D SETS UP A TWO DIMENSIONAL PLOT AREA OF THE 
Sete iplED SIZE IN INCHES.  CAEL SCMPLX SELECTIS THE LET- 
TERING TYPE FOR THE HEADING. CALL MESSAG AND CALL REALNO 
ALLOW ME TO PLACE A MESSAGE IN THE PLOT AREA AT A SPECIFIC 
EHI SPECIFIED IN INCHES IN THE CALL. CALL  XNAME AND 
YNAME ALLOW ME TO LABEL THE AXES. CALL HEADIN ALLOWS ME 
TO PLACE A TITLE ON THE PLOT. CALL THKFRM AND FRAME ALLOW 
ME TO PUT FRAMES AROUND THE PLOT AREA AND THE WHOLE PLOT 
ECET THICKNESS. CALL GRAF SETS UP THE PHYSICAL ORIGINS 
AND MAX VALUES AND THE INCREMENTS FOR EACH AXIS. CALL 
MUEVE TELLS THE MACHINE TO PLOT THE DESIRED ARRAYS OF DATA 
ESINTS. CERES NEEDS EHNPSEDRHATESPHESLEIC PROT-AND GIVES A 
MEW PAGE FOR THE NEAT PLOT. 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
MALE 
CALL 
CALL 


Hei e 

BLOW RA iS) 

YAXANG (0) 

IES >) 

AREAZD GS. , ©.) 

EM PIA 

UScASUMMASSSP ,100,0.1,5:5) 
REALNO(M,-2, 'ABUT' , 'ABUT' ) 

HES ACO KGP 100 ABUT , ABUT ) 
HES ACG A NVOCUME=$ .10050.1.5.) 


CALL 
CALL 


REALNO(VOL,-2,'ABUT','ABUT') 
DE AGRM’ TOQC ABUT'.'ABUT') 


SIS 


CALCE 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


MATES 


CALL 
CALE 
CALL 
CALE 
CALL 
CALL 
CALE 
CALL 
CAEL 
CATE 
CALL 
CALL 
CALL 
CALL 
CAEL 


THUS 


CALL 
CALL 
CALL 
CALL 
CALL 
CALLE 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


THS 


ANAME( FREQUENCY (IN) HERTZ) teen 

YNAME ("VOLTAGE (IN VOLTS)$",100) 

HEADIN( "VOLTAGE VS. FREQUENCYS LO NM 
THK PERM (TORE 

FRAME 

GRAF (0.0, INC,MAX,7.5,YINC4,YMAX4) 

CURVE CDTI TOOTO 

ENDPC COJ 


SECTION PLOTS PISTON VELOCITY VERSUSOENEOUEDNSM 


AREAZO Seer) 

MESSAG('MASS-$',100,0.1,5.5) 

REALNO(M,-2, ABUT', 'ABUT ') 
MESSAG('KG$',100,'ABUT', ABUT') 

MESSAGE VOLUME=SS ADO OA 
REALNO(VOL,-2. ABUT’ , ‘ABUT’ ) 
MESSAG('M**3$S' TOO “ABUT | Seu 

ANAME( 'FREQUENCY(IN HERTZI)S". LON 

YNAME( 'PISTON VELOCITY(IN METERS/ SECS ieee 
HEADIN('PISTON VELOCITY VS. FREQUENCY$ TOO MB 
THRP RMO?) 

FRAME 

GRAP (O OVINO MA COTO TINC EAA 
CUEVE CDAS = 00.105 

ENE ELECO 


SECTION PLOTS ACOUSTIC POWER VERSUS FREGUE IG 


AREAZD(8..6.) 
MESSAGE MASS => OOO SON 
REALNO(M{—-2> ABUL 2 ABUT) 

MESSAG( KGS’ ,100, ABUT = Asa) 

MESSAG( VOLUME=3" TULA 
REALNO(VOL,-2, 'ABUT’ ,'ABUT') 

MESSAG( M*"Sp' 100 ABUT T ABUTIS 
XNAME(' FEREQUENCY CIN HERTAZ)S COUP 
YNAME( ' ACOUSTIC POWER (IN N-M/S)$',100) 
HEADIN('ACOUSTIC POWER VS. FREQUENCYS 01000 M 
TANIA SS 

FRAME 

GRAF(O0O.O,INC,MAX,0-O-YINCS N 
CURVECDSEH. 3000109) 

ENDPC (0) 


SECTION PLOTS POWER RATIO VERSUSSSEEGQUEDLSIS 
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CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


CALL 


CALL 
STOP 
END 


AREAS O.) 

PES C A V Uno ol O0,0,1,5.5) 
REALNO(M,-2,'ABUT' , 'ABUT' ) 

PESSAC Nee 100. ABUT, “ABUT ) 
Ve Sono VOLUMES 100,0.1,5.) 
REALNO(VOL,-2,'ABUT', 'ABUT ') 

MES ACI SOS ABUT-, ABUT ' ) 
XNAME('FREQUENCY(IN HERTZ)$',100) 


ENARE T ACOUSTIC PWR/ ELECTRIC PWRE .100) 
HEADIN(' POWER RATIO VS. FREQUENCY$' ,100,2. 


THKERM( . 02) 

FRAME 
GRAPRMOROTTITNC MAX. 0.0, YINC3,YMAX3) 
CURVE (D O73000) 

ENDPL(0) 


TL) 


DONEPL ENDS THE PLOTTING SEQUENCE AND SENDS YOU 


Back TO CMs. 


DONEPL 


>> 


APPENDIX C. STAR DRIVER DIMENSIONAL DRAWINGS 





Figure C-1. Assembly Drawing of NIB Driver with 
all Housing Parts and Test Lid. 


156 





DRIVER HOUSING - RESONATOR VIEW 


(Dimensions in inches ) 


Figure C-2. Driver Housing--Resonator View. 
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Figure C-3. Driver Housing--Side View Section AA 


Showing Electrical Leads, Bolt Holes, 
Microphone Penetration, and FET 
Follower Housing. 


158 


| 


and FET location. 





VIEW AA 


AS IST e C» 


(Dimensions in inchee ) 


Figure C-4. Feed-Through for Electrical Leads 
Shown in Figure V-5. 
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Figure C-5. Driver Housingz--Side View Section BB 
Showing Pressure Gage and Helium 


Fill Ports3 
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4 Equally spaced 
clearance holes for 


19/32 bolts at 


Te 





VIEW AA 


eC Eo OURE TRANSDUCER 


(Dimensions in inches) 


Figure C-6. DC Pressure Gage Feed- 
Through Plug Details. 
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Figure C-7. 





DRIVER HOUSING - LIB iS 


Driver Housing--Lid View Showing 
Capillary, Microphone, and Follower, 
and Bolt Circles for Resonator 
Attachment, Pusher Plate, Pressure 
Lid and GAS Can Lid. 





(Note: "Plete is 3/18 wide end 


end slopes l0" et outer 


edge. ? 


HS E E TO VIEW) 


(Dimensions in inches). 


Figure C-8. Aluminum Pusher Plate. 
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Groove for 


1/16" lesd 
O-ring e 


Y 


(Note: O-ring groove U.018 wide 
and .@36° deep: 


lid 1s 1/2” wide] 


PRESSURE LID IDA 


(Dimensions in inches) 


Figure C-9. Pressure Lid with Lead 
O-ring Groove. 
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Pressure and Helium Leak Test Lid. 


Figure C-10. 
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B Enually eoeced 
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4 Equelly epeced 
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REDUCER CONE 


(Dimensione in inches ) 


Figure C-11. Reducer Cone. 


166 


ARREND AP TAR DRIVER LIST ORE MATERIALS 


The driver consists of a modified JBL model 2450J neo- 
dymium-iron-boron magnet with iron pole pieces on either 
side. The diaphragm is titanium and the voice coil is 
aluminum. 

The housing, pusher plate, pressure lid, and reducer 
cone will all be machined from 6061-T6 aluminum. All open- 
Peo the dager housing will be sealed with 1/16" or 
1/32" diameter lead o-rings. The electrical  feed-throughs 


will be sealed with Emerson & Cuming's STYCAST «tm» 2acoPr 


E 


epoxy. All wire leads are copper. The steel bolts connectin 


m 


the resonator to the driver are size 1/4-20. All other steel 


bolts are size 10-232. 


Doy 
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